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The theory of matrix models is reviewed from the point of view of its relation to integrable hierarchies. Determinantal formulas, relation to conformal 
field models and the theory of Generalized Kontsevich model are discussed in some detail. Attention is also paid to the group-theoretical interpretation of 
T-functions which allows to go beyond the restricted set of the (multicomponent) KP and Toda integrable hierarchies. 
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1 Introduction 

The purpose of these notes is to review one of the branches of modern string 
theory: the theory of matrix models. We put emphasize on their intrinsic 
integrable structure and almost ignore direct physical applications which are 
broadly discussed in the literature. Also most of technical details and references 
are omited: they can be found in a recent review [1]. 

Both "matrix models" and "integrability" are somewhat misleading names 
for the field to be discussed, they refer more to the history of the subject than 
to its real content. In fact the problem which is actually addressed is that 
of description of non-perturbative partition functions in quantum theory. The 
term "non-perturbative partition function" is now widely used to denote the 
generating functional of all the exact correlation functions in a given quantum 
model. Such quantity is given by a functional integral where the weight in the 
sum over trajectories is defined by effective action, which contains either all 
the possible (local or non-local) counterterms or generic coupling_Jjo external 
fields, so that any correlator can be obtained as derivative with respect to 
appropriate coupling constants or background fields. These exact generating 
functionals possess new peculiar properties, resulting from the possibility to 
perform arbitrary change of integration variables in the functional integral. Such 
properties are never studied in the orthodox quantum field theory because there 
the freedom to change integration variables is severely restricted by requirements 
of locality and renormalizability, which at last lost their role as fundamental 
principles of physics with creation of string theory. 

Every change of integration variables can be alternatively described as some 
change of parameters of non-perturbative partition function (i.e. the coupling 
constants or background fields in effective action). Thus invariance of the in- 



tegral implies certain relations (Ward identities) between partition functions at 
different values of parameters. Since all the fields of the model are integrated 
over, the set of relations is actually exhaustively large: more or less any two sets 
of parameters are related. Exact formulation of this property is yet unknown. 
The natural first step in these investigations is to look at the finite-dimensional 
integrals and then proceed to functional integrals by increasing the number of 
integrations. In turn, the natural way to do it is to make use of the well studied 
"matrix models" , which deal with N x N matrix integrals and their behaviour 
in the large- N limit. It appears that non-perturbative partition functions of 
matrix models, at least when they can be handled with the presently avail- 
able techniques, are closely related to "r-functions" , introduced originally in 
the study of integrable hierarchies. It looks very probable that some crucial 
characteristics of such partition functions are in fact not so peculiar for these 
simple examples, but remain true in the absolutely general setting. Extraction 
of such properties and construction of the adequate notion of "generalized r- 
function" is the main task of further work in the theory of "matrix models" and 
"integrable systems" . 

One of the most straightforward and still promissing approaches is based 
on interpretation of Ward identites for non-perturbative partition functions as 
Hirota-like equations (supplemented by a much smaller set of "string equa- 
tions"), while generalized r-functions, which are solutions to these equations, 
are interpreted as group-theoretical objects (generating functionals of all the 
matrix elements of a group clement in particular representation). If successfull, 
this approach can provide description of exact correlation functions in terms 
of some (originally hidden) symmetry of the given class of theories, thus rais- 
ing to the new height the relation between physical theories and symmetries, 
which was the guiding line for developement of theoretical physics during the 
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last decades. 

Analysis of non-perturbative partition functions is very important for one 
more reason. Construction of a generating functional is essentially "exponen- 
tiation of perturbations" , i.e. it deformes original (bare) action of the model. 
When perturbation parameters (extra coupling constants or background fields) 
are non-infinitesimal, one in fact obtains entire set of models instead of origi- 
nal one. Moreover, original model is no longer distinguished within this cl asa : 
non-perturbative partition functions are associated with classes of models, feet 
with a single model. One can easily recognize here realization of the main idea 
of the string programm (see for example [2]). The study of non-perturbative 
partition functions even for such a simple class as ordinary matrix models can 
lead to much better understanding of the general idea. As usual, this can help 
to figure out what the adequate questions are and to develop effective technique 
to answer these questions. 

The purpose of these notes is to briefly illustrate the general ideas with very 
simple examples. A lot of work is still required to obtain applications to the 
really interesting problems. Still, simple examples are enough to understand 
the ideas, and often conceptual level is no less important than that of technical 
effectiveness. 

We begin in the following sections from cosideration of the simplest matrix 
models, to be refered to as discrete and Kontsevich (continuous) models. Again, 
these names reflect more the history than the real content of the subject (con- 
tinuous models were originaly described non-explicitly as specific (multiscaling) 
large- N limits of the discrete ones). In the context of non-perturbative partition 
functions the difference is that discrete models possess effective actions with all 
the possible counterterms added, while in Kontsevich model an external field 



(source) is introduced. Analysis of these theories includes their characterization 
as eigenvalue models and derivation of related determinant representations. We 
also consider description of discrete models in the language of conformal field 
theory. It is important for connecting these matrix models to the physically 
relevant Liouville theory of 2c? gravity and - more essential for our notes - to 
the concept of KP and Toda r-functions. 

Then we turn to the free-fermion (Grassmannian) description of KP/Toda r- 
functions and list some results from the theory of Generalized Kontsevich model 
[3]-[8]. Some other interesting models, which look a priori very different, are 
in fact just particular examples of Kontsevich model, which actually describes 
a big family of theories. This example should be very instructive for the fu- 
ture understanding of the interplay between perturbative and non-perturbative 
information contained in the non-perturbative partition function. 

The last topic of these notes concerns group-theoretical interpretation of r- 
functions. The natural object, which arises in this way, while possessing the 
most important properties of conventional r-functions, is in fact much more 
general. In this framework the KP/Toda r- functions are associated with funda- 
mental representations of SL{N) and the closely related theory of the simply- 
laced Kac- Moody algebras of level k = 1. In fact r- function can be easily defined 
for any representation of any group (including also quantum groups - this can 
be important for the future construction of string field theory, where the idea 
of "third quantization" requires consideration of operator-valued r-functions; 
since it takes our non-perturbative partition function as input - effective action 
of the full string field theory). 



2 The basic example: discrete 1-matrix model 



The sample example of matrix model is that of 1-matrix integral 

Z N {t} = c N [ dHe^7=o tkTlHk , (2.1) 

JNxN 

where the integral is over N x N matrix H and dH = J\ i j dHj. 1 This measure 
is invariant under the conjugation H — > U HIP with any unitary N x N matrix 
U, and the "action" J2T=o tk^cH in (2.1) is the most general one consistent 
with this invariance. Thus Z^{t} is indeed an example of the non-perturbative 
partition function in the sence, described in the Introduction. All observables 
in the theory are given by algebraic combinations of TrH k and their correlation 
functions can be obtained by action of i^-derivatives on Zpj{t}. Our goal now 
is to hnd some more invariant description of this quantity, not sohspecific as the 
matrix integral (2.1). 

2.1 Ward identities 

Such description is provided by Ward identities. The integral is invariant ] ur[ dcr 
any change of integration matrix- variable H — > f(H). It is convenient to choose 
the special basis in the space of such transformations: 

5H = e n H n+1 . (2.2) 

1 This integral is often refered to as Hermitean. In most of our considerations we do not 
need to specify integration contours in matrix integrals, in particular eigenvalues of Hij do not 
need to be real. What is indeed important is the "flatness" of the measure dH = Y\ i j dHij. 
Below "Hermitean" (as opposed, for example, to "unitary") will imply just this choice of the 
measure and not any reality condition. 



□ 



Here e„ is some infinitesimal matiix and, of course, n > — 1. Invariance of the 



integral implies the following identity: 

/ rf// e E:„^ k = /^ff+eX+^SL^^^ 1 )^ (2 .3) 

JNxN J 

i.e. 



= 0. (2.4) 

In order to evaluate the Jacobian Tr SH SI ^ let us restore the matrix indices: 

n n 

{5H n+1 ) l3 = Y,(H k SHH n -% = J2(H k h(5H) lm (H n ~ k ) mj , (2.5) 

k=0 k=0 

and to obtain Tr ^jjj- put I = i and m = j, so that 

Tr SH = TrH k TrH n ~ k . (2.6) 

fc=0 

Any correlation function can be obtained as variation of the coupling constants: 
< TrH ai ...TrH a " >= J dHe^T=o tkTtHk TrH ai ...TrH a *> = 



Qn 

-Z N {t}. 



(2.7) 



dt ai ...dt a „ 

This relation together with (2.6) can be used to rewrite (2.4) as: 

L n Z N {t} - 0, n > -1 (2.8) 

with 

°° d " d 2 
L n = Y,kt k - +V-— . (2.9) 

Note that according to the definition (2.1) 

^-Z N = NZ N . (2.10) 



2.1.1 Details and comments 



Several remarks are now in order. 

First of all, expression in brackets in (2.4) represents just all the equations of 
motion for the model (2.1), and (2.8) is nothing but another way to represent 
the same set of equations. This example illustrates what "exhaustively large" 
set of Ward identities is: it should be essentially the same as the set of all 
equations of motion. 

Second, commutator of any two operators L n apearing in (2.8) should also 
annihilate Z/v{£}. Another indication that we already got a complete set of 
constraints, is that L n 's form a closed algebra: 

[L n ,L m ] = (n-m)L n+m , n, m>-\. ( 2 -ll) 

□ 

Its particular representation (2.9) is refered to as "discrete Virasoro algebra" (to 
emphasize the difference with "continuous Virasoro" constraints, see eq.(3.64) 
below). 

Third, (2.8) can be considered as invariant formulation of what is Zn'- it is a 
solution of this set of compatible differential equation s. F rom this point of view 
eq.(2.1) is rather a particular representation of a t is sensible to look for 
other representations as well (we shall later discuss two <j)TT|licm: one in terms 
of CFT, another in terms of Kontsevich integrals). 

Fourth, one can try to analyze the uniqueness of the solutions to (2.8). If there 
are not too many solutions, the set of constraints can be considered complete. 
A natural approach to classification of solutions to the algebra of constraints 
is in terms of the orbits of the corresponding group [9]. Let us consider an 
oversimplified example, which can still be usefull to understand implications of 
the complete set of WI as well as to clarify the meaning of classes of universality 



and of integrability. 

Imagine, that instead of (2.8) with L„'s defined in (2.9) we would obtain 
somewhat simpler equations: 2 

l n Z = 0, n > with l n = kt k - . (2.12) 

k=i dtk+n 

Then operator li can be interpreted as generating the shifts 

*2 — ► *2 + £1*1) 

t 3 ^t 3 + 2e 1 t 2 , (2.13) 



□ 



We can use it to shift t 2 to zero, and eq. Z|=^p then implies that 

Z(h,t 2 ,t 3 ,...) = Z(tx,0,t 3 ,...) (2.14) 

(h = t k - k > 3). 1=1 

Next, operator l 2 generates the shifts 

*3 — * *3 + £2*1, 

U — >t 4 + 2e 2 t 2 , (2.15) 

□ 

and does not affect t 2 . We can now use eq. l 2 Z — to argue that 

Z{t 1 ,t 2 ,t s ,U,...)=Z(t u 0j 3 ,U,...) = Z(t 1 ,0,Q~U,---) ( 2 - 16 ) 



2 One can call them "classical" approximation to (2.8), since they would arise if the variation 
of measure (i.e. a "quantui p ep ect") was not taken into account in the derivation of (2.8). 
Though this concept is often used in physics it does not have much sense in the present 



context, when we are analyzing exact properties of functional (matrix) integrals 



etc. Assuming that Z is not very much dependent on tk with k — > oo, 3 we can 
conclude, that 

Z(h,t 2 ,t 3 ,...) = Z(h,0,0,...) = Z(1,0,0,...) (2.17) 

(at the last step we also used the equation IqZ = |Q — tj o rescale t\ to unity). 

All this reasoning had sense provided t\ 7^ 0. Otherwise we would get 
Z(0, 1,0,0,...), if t x = 0, t 2 £ 0, or Z(0, 0,1,0,...), if t x = h = 0, i 3 ^ 
etc. In other words, we obtain classes of universality (such that the value of 
partition function is just the same in the whole class), which in this oversimpli- 
fied example are labeled just by the first non-vanishing time- variable. Analysis 
of the orbit structure for the actually impgf^ant realizations of groups, like that 
connected to eq.(2.9) has never been performed in the context of matrix model 
theory. 

In this oversimplified case the constraints actually allow one to eliminate all 
the dependence on the time-variables, i.e. to solve equations for Z exactly. In 
realistic examples one deals with less trivial representations of the constraint 
algebra, like (2.11). It appears that in this general framework constraints some- 
how imply the integrability structure of the model, what can thus be considered 
as a slightly more complicated version of the same solvability phenomenon. 

2.2 CFT interpretation of 1-matrix model 

Given a complete set of the constraints on partition function of infinitely many 
variables which form some closed algebra we can now ask an inverse question: 
how these equations can be solved or what is the integral representation of 

3 This, by the way, is hardly correct in this particular example, when the group has no 
compact orbits. 



partition function. One approach to this problem is analysis of orbits, briefly 
mentioned at the end of the previous subsection. Now we turn to another 
technique [10], which makes use of the knowledge from conformal field theory. 
This constructions can have some meaning from the "physical" point of view, 
which implies certain duality between the 2-dimensional world surfaces and 
the spectral surfaces, associated to configuration space of the string theory. 
However, our present goal is more formal than discussion of this duality: we are 
going to use the methods of CFT for solving the constraint equations. 

This is especialy natural when the algebra of constraints is Virasoro algebra, 
as is the case with the 1-matrix model, or some other algebra known to arise as 
a chiral algebra in some simple conformal models. In fact the approach to be 
discussed is rather general and can be applied to construction of matrix models, 
associated with many different algebraic structures: the only requirement is 
existence of the (massless) free-field representation. 

We begin from the set of "discrete Virasoro constraints" (2.8). The CFT 
formulation of interest should provide the solution to these equations in the 
form of some correlation function in some conformal field theory. Of course, 
it becomes natural if we somehow identify the operators L n (2.9) with the 
harmonics of a stress-tensor T n , which satisfy the same algebra, and manage 
to relate the constraint that L n annihilate the correlator to the statement that 
T„ annihilate the vacuum state. Thus the procedure is naturally split into two 
steps. First, we should find a i-dependent operator (" Hamiltonian" ) H(t), such 
that 




(2.18) 



This will relate differential operators L n to T n 's expressed through the fields of 
conformal model. Second, we need to enumerate the states, that are annihilated 



5 



by the operators T„ with n > — 1, i.e. solve equation 



It is easy to check that 



T n | G) = 



(2.19) 



□ 



for the ket-states, what is an internal problem of conformal field theory. If both 
ingredients H(t) and | G) are found, solution to the problem is given by 



G). 



(2.20) 



To be more explicit, for the case of the discrete Virasoro constraints we can 
just look for solutions in terms of the simplest possible conformal model: that 
of a one holomorphic scalar field 



(z) = q + p\ogz + 



J-k k 



(2.21) 



[Jn, Jm] = nS n+m ,0, [q,P] = 1- 

Then the procedure is as follows. Define vacuum states 

J fc |0)=0, (iV|J_ fe = 0, fc>0 
p|0) = 0, (N\p = N(N\, 

the stress-tensor 

T{z) = \[d^{z)] 2 =^T„z-"- 2 , T n = Y,J-hJk+n + \ JaJb 



(2.22) 



fc>0 



a,6>0 



(2.23) 



and the Hamiltonian 



ff(*) = 4=X>.7* = 4= / U(z)J(z) 
^ 2 fc>o V2 Jc ° 

C/(z) = ^i fe z fe , J(z) = d<l>(z). 

k>0 



(2.24) 



and 



L n (N\e H M = (N\e H ^T n 

□ □ 



T„|0) = 0, n > -1. 
As an immediate consequence, any correlator of the form 

Zjv{t | G} = (JVle^WGlO) 

gives a solution to (2.8) provided 

[T n ,G] = 0, n>-l. 



(2.25) 



(2.26) 



(2.27) 



(2.28) 



In fact operators G that commute with the stress tensor are well known: these 
are just any functions of the "screening charges" 4 



J± 



(2.29) 



The correlator (2.27) will be non-vanishing only if the matching condition for 
zero-modes of <p is satisfied. If we demand the operator to depend only on 
Q_l_, this implies that only one term of the expansion in powers of Q + will 
contribute to (2.27), so that the result is essentially independent on the choice 
of the function G(Q+), we can for example take G(Q+) — and obtain: 

1 



Z N {t}~—(N\e H W(Q + ) N \0). 



(2.30) 

This correlator is easy to evaluate using Wick theorem and the propagator 



4 For notational simplicity we omit the normal ordering signs, in fact the relevant operators 
are : e H : and : : 
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(z)4>(z') ~ log(z — z'). Finally we get 



7 U\ 1 I AT I 75 £ U(z)d<j>(z) 



»=i 



o) 



□ 

(2.31) 



in the form of a multiple integral. This integral does not yet look like the matrix 
integral (2.1). However, it is the same: (2.31) is an "eigenvalue representation" r— | 
of matrix integral, see [11] and eq.(2.53) in the next subsection 2.3. 

2.2.1 Details and comments 

Thus in the simplest case we resolved the inverse problem: reconstructed an 
integral representation from the set of discrete Virasoro constraints. HoOer, 
the answer we got seems a little more general than (reflmamo) and (2.53): the 
r.h.s. of eq.(2.31) still depends on the contours of integration. Moreover, we 
can also recall that the operator G above could depend not only on Q +1 but 
also on Q_. The most general formula is a little more complicated than (2.31): 

' -(N\e H ^(Q + ) N + M (Q.) M \0) = □ 



Z N {t | Cj, C r } 



1 



(N + M)\M\ 



(N + M)\M\ 

N+M , M „ 

(2.32) 



N+M M 

TT J> d Zi e u ^ TT i dz' r i 

i=l J Ci r=l J C'r 



t-iN+M, 
ili<i v 



i<j \^ ) LLr<8 



nf +M nf(^-^) 2 



□ 



We refer to the papers [10] for discussion of the issue of contour-dependence. In 
certain sense all these different integrals can be considered as branches of the 
same analytical function Z^{t}. Dependence on M is essentially eli riink teH by 
Cauchy integration around the poles in denominator in (2.32). 



Above construction can be straightforwa rdly applied to any other algebras of 
constraints, provided: 

(i) The free-field representation of the algebra is known in the CFT- 
framework, such that the generators are polinomials in the fields <f> (only in 
such case it is straightforward to constru Hamiltonian H, which relates 
CFT-realization of the algebra to that in terms of differential operators w.r.to 
the t- variables; in fact under this condition H is usually linear in £'s and </>'s). 
There are examples (like Frenkel-Kac representation of level k = 1 simply-laced 
Kac- Moody algebras [12] or generic reductions of the WZNW model [13], [14]- 
[17]) when generators arb exponents of free fields, then this construction should 
be slightly rhodtified. 

(ii) It is easy to find vacuum, annihilated by the relevant generators (here, for 
example, is the problem with application of this approach to the case of "con- 



tinuous" Virasoro and VF-constraints). The resolution to this pro )lexi involves 



consideration of correlates on Riemann surfaces with non-trivial topologies, of- 
ten - of infinite genus. 

(iii) The free-field representation of the "screening charges", i.e. operators 
that commute with the generators of the group within the conformal model, is 
explicitly known. 

These conditions are fulfilled in many case in CFT, including conventional 
W-algebras [18] and AT — 1 5 supersymmetric models. 

For illustrative purposes we present here several forrrrmas from the last paper 



5 In the case of . V 



screening charges. 



7tr 



: 2 supcrsymmctry a problem arises because of the lack of reasonable 
the most naive level the relevant operator to be integrated over su- 



perspace (over dzd^O) in order to produce screening charge has dimension 1 — ^M, which 
vanishes when M = 2. 
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of ref. [10] for the case of the W r+ i-constraints, associated with the simply-laced 
algebras Q of rank r. 

Partition function in such "conformal multimatrix model" is a function of 
" time- variables" vt , k — O...00, A = l...r = rankC/, and also depends on 
the integer-valued r-vector TV = {Ni...n r }. The W r+ i-constraints imposed on 
partition function are: 

Wi a) (t)Z^{t} = 0, n>l-a, o = 2...r + l. (2.33) 

The form of the M^-operators is somewhat complicated, for example, in the case 
of r + 1 = 3 (i.e. for Q = SX(3)) 



j-f K dt a dt h ^ dt a dt b ' 

a+o—n 



(2.34) 



W n } = Yl WU-aA kt k lti^A 2kt k lii-J- )+ 

£jT Q Otk+n+l Otk+n+l Otk+n+l 



fc>o 



V- n 92 , 92 n," d2 



dt a dt b dt a dt b ' dtadh) 

.a-\-b—n-\-K 



(2.35) 



4 >p d 3 d 3 



3 , 77* dt a dt b dt c dt a dt b dt c 

a+b+c—n 

and two types of time- variables, denoted through t k and tu- are associated with 
two orthogonal directions in the Cartan plane of A2 : e = ^= , e = v ^| 2 . 6 

All other formulas, however, are very simple: Conformal model is usually that 
of the r free fields, S ~ J 8<j)d(j)d 2 z, which is used to describe representation of 

6 Such orthogonal basis is especially convenient for discussion of integrability properties of 
the model, these t and t are linear combinations of t ime- variables i» appearing in eqs. (2.36) 
and (2.41). ' ' 



the level one Kac- Moody algebra, associated with Q. Hamiltonian 

ff(^..^ +1) H££4 A W fe , (2.36) 

a=i k>a 

where are associated with "fundamental weight" vectors v\ in Cartan 

hyperplane and in the simplest case of Q = SL(r + 1) satisfy 

j r+l 

P>\ ■ #v = S\\> ~~r, fix = 0, 

A=l 

thus only r of the time variables t^ . . . £( r+1 ) are linearly independent. Relation 
between differential operators W^it) and operators in the CFT is now 

defined by 



Wi a \N\e H ^ =(N\e H ^W ( i a \ 
a = 2, . . . ,p; i > 1 — a, 



(2.37) 



where 



W n a) = j> z a+ "- 1 W (a )(2) 



(2.38) 



A 

are spin-a generators of the W^ +1 algebra. The screening charges, that com- 
mute with all the (z) are given by 

Q(«) = I j( Q ) = /e 5,? (2.39) 



{a} being roots of finite-dimensional simply laced Lie algebra Q . 
Thus partition function arises in the form: 



Z<jL{t} = (N\e H ^G{Q^}\0) (2.40) 



where G is an exponential function of screening charges. Evaluation of the 
free-feild correlator gives: 



m - J n 



i=l \A;fc>0 

N a N/J 



(2.41) 



nniK* 



(«) _ z (/3))S/3 



(a,/3) i=lj = l 

In fact this expression can be rewritten in terms of an r-matrix integral - a 
" conformal multimatrix model" : 



JNxN 1 



Yl Det (# (q) O J - I ® 



50 



(2.42) 



(«,/3) 



In the simplest case of W3 algebra eq.(2.41) with insertion of only two (of the 
six) screenings Q ai and Q a2 turns into 



Z NuN 2 (*) *) 



1 



iVi!AT 2 ! 



(iVi , iV 2 1 e ff (*■*) (Q( Q1 ) ) Wl (Q( Q2 ) ) N2 10) 



1 



N^.NJ. 



n / dZie^JJ / dy J e ^A(x)A(x,y)A(y), 



(2.43) 



where A(x,j/) = A(x)A(y)Y[ i j(xi — yj). This model is associated with the 
algebra Q = SL(3), while the original 1-matrix model (2.30)-(2.32) - with Q = 
SL(2). 

The whole series of models (2.41-2.42) for Q = SL(r + 1) is distinguished 
by its relation to the level k = 1 simply-laced Kac- Moody algebras. In this 



particular situation the underlying conformal model has integer central charge 
c = r = rank Q and can be " fermionized" . The main feature of this formulation 
is that the Kac-Moody currents (which after integration turn into "screening 
charges" in the above construction) are quadratic in fermionic fields, while they 
are represented by exponents in the free-boson formulation. 

In fact fermionic (spinor) model naturally possesses GL(r + 1) rather than 
SL(r + 1) symmetry (other simply-laced algebras can be embedded into larger 
Gi-algebras and this provides fermionic descriprion for them idie case of 
k = 1). The model contains r + 1 spin-1/2 fields tpi and their conjugate tpi 
(b, c-systems) ; 

r+l - 

S = J2 -4>]^]d 2 z, (2.44) 
3=1 

central charge c = r + 1 , and operator algebra is 

6 



z — z 



^{z)Mz') = (z - z')5 jk : ipj(z)Mz') : + (1 - S jk ) : ^(z)M^') ■ (2 ' 45) 
ipj(z)ipk(z') = {z - z')S jk ■ ^j{z)^ k {z') : + (1 - 5jk) ■ ^}j{z)^} k (z') : 



The Kac-Moody cur 



ents of the level-one GL(r + l) fc=1 are j ust Jjfe =: ipjipk '■ 



j, k = 1 . . . r + 1, and screening charges are — iE^' § : ip 



where E- k 
of GL{r + 1). 



are representatives of the roots a in the matrix representati : 
Cartan subalgebra is represented by Jjj , while positive and negative Borel sub- 
algebras - by Jjk with j < k and j > k respectively. In eq.(2.32) Q + = 

7 This is possible only for very special Kac-Moody algebras, and such formulation is im- 
portant in order to deal with conventional formulation of integrability, which usually involves 
commuting Hamiltonian flows (not just a closed algebra of flows) and fermionic realization 
of the universal module space (universal Grassmannian) . In fact these restrictions are quite 
arbitrary and can be removed (though this is not yet done in full details), see [1] and sections 
4, 5 below for more detailed discussion. 



i § i>i^2, Q-=i§ ifaipi while in eq.(2.43) Q {ai) = i § ipiifa, Q {oi2) = i § ^1^3 
(and Q<°») = iffafo, Q (ai) = i§^x, Q (qs) = i/^i, <2 (ae) = 
i / 413^2)- Q^" 6 ^ can be substituted instead of Q(" 2 ' in (2.43) without changing 
the answer. For generic r the similar choice of "adjacent" (not simple!) roots 
(such that their scalar products are +1 or 0) leads to selection of the following 
r screening operators — i § ^1^2 = —i § ^2^3, = i § ^3^4j • ■ -i 

i.e. Q«) = i § ?pjipj + i for odd j and Q^' — —i § ijjjipj+i for even j. 



2.3 1-matrix model in eigenvalue representation 

In the last section we found solution of Virasoro constraints (2.8) in the form 
of the multiple integral (2.31). Now we shall see that this expression is in 
fact equal to original matrix integral (2.1) and arises after "auxiliary" angular 
variables are explicitly integrated out. These angular variables are in fact the 
ones associated with physical vector bosons andl tae possibility to solve this 
sector of the model explicitly is peculiar feature of the simplest class of matrix 
models, naturaly named eigenvalue models. The theory of eigenvalue models 
is in a sense equivalent to the theory of conventional intcgrable hierarchies 



and thus is rather straightforward to work e n. Inclusion of non-trivial angular 
integration in the general scheme is still a sophisticated task, with no universal 
solution found so far. Also unknown is solution of inverse problem: how can be 
arbitrary eigenvalue model - with integration over eigenvalues of some matrix 
- "lifted to" the full matrix model where integration goes over entire matrix 8 ; 
in other words what is the way to couple vector bosons to the "topological" 
eigenvalue sector so that the two sectors are interacting only in the "solvable" 
fashion and angular integrations can be easily performed. 



s see also ref. [19] 



Let us now turn to particular example of the 1-matrix integral (2.1). First 
of all, this model possesses gauge symmetry, associated with the unitary (an- 
gular) rotation of matrices, H a — > U^H a U a . This illustrates the general phe- 
nomenon: matrix models are usually gauge theories. In the case of eigenvalue 
models this symmetry is realized without "gauge fields" V a p, which would de- 
pend on pairs of indices a, (3 and transform like V a p — > U^VapUp. In other 
words, eigenvalue models are gauge theories without gauge fields, i.e. are pure 
topological. The case of the 1-matrix model (2.1) 



Z N {t} = c N [ dHe^=o tkTlH \ (2.46) 

JNxN 



is especially simple, because separation of eigenvalue and angular variables does 
not involve any information about unitary-matrix integrals. Take 

□ 



□ 



H = U^DU, (2.47) 

□ 

where U is a unitary matrix and diagonal matrix D = diag(/ii.../i7v) has eigen- 



10 



values of H as its entries. Then integration measure 9 



□ 



N r , T n N 



dH=f[ dH l3 = \ dU] f[dhiA 2 (h), (2.51) 

i jtl [dUcartanl f_\ 

where " Van-der-Monde determinant" A(h) = det^h{~ 1 — Y[i>j(^i ~ hj) an d 
[dll] is Haar measure of integration over unitary matrices. 

It remains to note that the "action" Tr U(H) = Y^k=a^k^H k with H sub- 
stituted in the form (2.47) is independent of U; 

N 

Tr U(H) = ^2u(h l ). (2.52) 

i=l 



9 In order to derive eq. (2.51) one can consider the norm of infinitesimal variation 

JV JV 

|| sh \\ 2 = ^2 I m a l 2 = SH v SH a = T < SH ) 2 = 

, (2-48) 
= Tr (-WSUIPDU + D8U + U 1 SDU) = 

= Tr(<5D) 2 + 2iTrSu[8D, D] + 2Tr (-SuDSuD + {8uf D 2 ) , 

where 8u = jSUU^ = 8vJ and 8D = diag(<5/ii . . . <5/ijv). The second term at the r.h.s. vanishes 
because both D and 8D are diagonal and commute. Therefore 

JV JV 

|| SH || 2 = ^(Sh,) 2 + ^ (5«)ij(<5«)ii(h« - h 3 ) 2 . (2.49) 

i — 1 — l 

Now it remains to recall the basic relation between the infinitesimal norm and the measure: 
if || 81 \\ 2 = G ab 8l a Sl b then [dl] = y/det ab G ab FJ dl a , to obtain eq.(2.51) with Haar measure 
[dU] = Y\ij d-Uij being associated with the infinitesimal norm ' ' 

JV JV 
|| Su \\ 2 = Tr(<5«) 2 = SuijSuji = | Suij \ 2 (2.50) 



i,j=l 

and [dUcartan] = Yl i=1 du ii ■ 



(2.53) 



Thus 

1 1 JV . JV 

z ^ = w\ n J dh ^ u{iu) n> " = 

i=X i>j 
i—l 

provided cjv in (2.1) and (2.46) is chosen to be 

CN=Nl W^W' (2 - 54) 

where the volume of unitary group in Haar measure is equal to 

Vol uW = ' . (2.55) 

ru=i k - 

2.4 Kontsevich-like representation of 1-matrix model 

Matrix-integral representation (2.1) is, howevejrTJtot the only possible one for 
the given eigenvalue model (2.31). Expression (2.1) involves the "most general 
action" , consistent with the symmetry H — ► UHU' . As was already mentioned 
in the introduction, alternative representation of the partition function should 
instead involve the general coupling to background (source) field. In the theory 
of matrix models such representations are known under the name of Kontsevich 
models. They will be the subject of detailed discussion in the next sections of 
this paper. What we need now is one simple identity, relating original (2.1) and 
Kontsevich-like representations of the 1-matrix theory: 

□ 



^jv{io = 0;i* = -^rA- k + ±4,2} J Nx N dH ^=° tk 



JJVxJV (25g ^ 

/ dX(detX) N e~ tl2 £ +tlAX = Zx2_{N,t}, 

Jnxn 2 



(27r)T-(detA) W Jnxn 



1 



□ 



where Z^{tk — 5^.2} = ( — 2-7r) ™ c/y. This relation follows from another iden- 
tity: 

J NxN dHe^ H2 Det{K ® I - 1 ® H) _ 

J ™ (2.57) 

L, n dXe-^ X2 det N (X + A) 



f v dX 

Jnxn 



e 2 



trX 2 



which is valid for any A and can be proved by different methods: see [1] and 
references therein. Note that integrals on the right and left hand sides are of 
differents sizes: NxN at the l.h.s. and nxn at the r.h.s. While TV-dependence is 
explicit at both sides of the equation, the n-dependence at the l.h.s. e ptep s only 
implicitly: through the allowed domain of variation of variables tk — — jtrA~ k + 
5<5/c.2- This is the usual feature of Kontsevich integrals: explicit n-depenence 
disappears once the integral is expressed through the t-like variables. 

Eq.(2.56) can be used to perform analytical continuation in N and define 
what is Zn for N, which are not positive integers. Since cn = for all negative 
integers (see [1]), the same is true for Z^. This property (tn — for all negative 
integers TV) is in fact characteristic for r-functions of forced hierarchies of which 
the partition function (2.53) is an example. 



3 Generalized Kontsevich Model (GKM) 

Let us now proceed to investigation of Kontsevich models of a rather general 
type. Further generalizations, leading directly to theories with physical vector 
bosons (see for example [20]), are beyond the scope of the present notes. The 
basic mathematical fact, responsible for solvability (integrability) of Kontsevich 
models is Duistermaat-Hcckmann theorem, which allows to evaluate explicitly 



the celebrated non-linear Harish-Chandra-Itzykson-Zuber integral over unitary 
matrices and thus transform matrix integral into an eigenvalue model. 



3.1 Kontsevich integral. The first step 

Kontsevich integral is defined as 

T v . n {L}= f dXe- trV( - x ^ Nt ^ x +^ LX . (3.1) 

J nxn 

In fact it depends only on eigenvalues of the matrix L. Indeed, substitute 
X = U 1 x D x Ux; L = U{d l U l in (3.1) and denote U = U X U{. Then 

= f[[ dx a x»e- v ^A 2 (x) I J dU] exp ( £ xj b \ U aS 

„_1 J J nxn [u-UCartanl \ „ fe=1 



(3.2) 

In order to proceed further one needs to evaluate the integral over unitary 
matrices, which appeared at the r.h.s. 



This integral can actually be represented in two different forms: 

In{X,L}= f 
J n 



[dU] 

tvXULU 1 



[dUcartan} 

3.3 

[dm c r.:^*JM 

nxn 

[dU c artan\ 

(the U's in the two integrals are related by transformation U — > UxUU^ 
and Haar measure is both left End right invariant). Form CIal (3.3) implies that 



I n {X, L} satisfies a set of simple equations [21]: 



d Nfc 



tr(^-J - trL k J I n {X, L} = 0, k > 0, 

71 ' fe \ (3 ' 4) 
tr f«rM - tTXk I in{x,L} = o, k>o, 



^ \dL tr 



which by themselves are not very restrictive. However, another formula, (3.3), 
implies that I n {X, L} in fact depends only on the eigenvalues of X and L, and 
for such I n {X, L} = I{x a , lb} eqs.(3.4) become very restrictive 10 and allow to 
determine I{x a , h} unambigously as a formal power seria in positive powers of 
x a and lb- The final answer is 



I n {X, L\ 

Normalization constant can be defined bjQking L = 0, when 



(3.6) 



Vol u{n) _ (2tt) 



I n {X, L = 0} = '7 ft = ^ ,, , (3.7) 



and using the fact that 

det ab fa(W '" ' 



A(/) 



II ^ J dBd^UO). (3.8) 



3.2 Itzykson-Zuber integral and Duistermaat-Heckmann 
theorem 

Eq.(3.6) is usually refered to as the Itzykson-Zuber formula [22]. In mathemat- 
ical literature it was earlier derived by Kharish- ( pr^ ndra [23], and in fact the 
integral (3.3) is the basic example of the coadjoint orbit integrals [24]- [26], which 
can be exactly evaluated with the help of the Duistermaat-Heckmann theorem 



10 When acting on /, which depends only on eigenvalues, matrix derivatives turn into: 



d_ 

8X tr ^-^ dx a 



dXn 



tr 3 r - V —/ 1 f— — ^ /• 

^tr ^ & x l X a - X b \ dx a dx b J 



etc. 



[27]- [30]. We now iterrupt our discussion of Kontsevich model for a brief illus- 
tration of this important phenomenon. Duistermaat-Heckmann theorem claims 
that under certain reLfedctive conditions (that dynamical flow is consistent with 
the action of a compact group) some integrals can be expressed in a simple 
way through extrema of the integrand. This almost sounds like a statement 
that quasiclassical approximation can be exact, with two correction that all 
the extrema - not only the deepest minimum of the action - should be taken 
into account, and that the quantum measure should be adjusted appropriately. 
When applicable, the theorem states that 

Q2 S \ -1/2 



WK SW - E (W e~ SW (3-9) 

tj>: dS/dd>=0 V 7 

The simplest example is given by the jia^egral J* [sin 6d0] e M cos 9 — - . We 

shall not dwell upon the rea^oas why DH theorem is true in the case of the 
Itzykson-Zuber integral (the basic Requirement: existence of the compact group 
action - that of unitary group - is obviously fulfilled in this case). Instead we just 
evaluate the r.h.s. of (3.9) provided the l.h.s. is given by J[dU] exp (trXULW). 
Then equation of motion for U looks like 

[X,ULrf]=0 (3.10) 

We assume that X and L are already diagonal matrices. Then (3.10) has 
an obvious solution U = I (identity matrix), but it is not unique. For 
generic diagonal X, L the most general solution is givem by arbitrary per- 
mutation matrix P: U — P. The "classical action" on such solution is equal 
to tiXULW = ^2 a x a lpf a \, while the preexponential factors provide Van-der- 
Monde determinants A(x)A(l) in denominator and the sign factor (— ) p . Since 

o 



^(-) P e^j^2x a lp {a) ^j =dete Xa 



we immediately obtain the IZ formula (3.6). Unfortunately the Duistermaat- 
Heckmann theory is not yet well developed and even if the vacuum average is 
exactly calculable, it does not provide immediate prescription for evaluation of 
correlators (or, what is essentially the same, corrections to DH formula - when 
it is not exactly true - are not yet described in a universal way). The very 
important general technique of exact evaluation of non — Gaussian unitary- 
matrix integrals is now doing its first steps (see [31]-[34]). 



3.3 Kontsevich integral. The second step 

Now we turn back to the eigenvalue representation of Kontsevich integral (3.1). 
Substitution of (3.6) into (3.2) gives: 

TvAL) = f[ / dx b e~ v ^A(x)±det ab e x ^ = 

(3.11) 



(2*) 



Atf) 



f[ f dx b x£e- v ^ +Xblb A(x), □ 

where we used antisymmetry of A(x) under permutations of x^s in order to 
change ^jdet ab e Xalb for e^"> Xblb underpffl|e sign of the x integration. 

We can now use the fact that A(x) = &et a bx1~ l in order to rewrite the r.h.s. 
of (3.11): 

TvAL} = {2^ deta % N{h \ (3-12) 

where 

ip a (l) ee ( dxx a - 1 e- v{x)+lx , a > 1. (3.13) 



These functions ip(l) satisfy a simple recurrent relation: 



d>f> a -i ( d 



a-1 



«-^-lsl 41 (3 ' 14) 



with 



$(Z) = ^(/) = / dxe- v{x)+lx 



(3.15) 



This completes the transformation of Kontsevich integral to the fibrin of eigen- 
value model. 



3.4 "Phases" of Kontsevich integral. GKM as the "quan- 
tum piece" of Ty{L\ in Kontsevich phase 

j — (p nc of the natural things to do in the study of functions, defined in the integral 
form, is to investigate various asimptotics when the arguments tend to various 
distinguished limits. In the case of Kontsevich integral the arguments are just 
eigenvalues l a of the matrix L, while their distinguished values are either zero or 
positions of singularities of potential V(x). For simplicity we assume that V(x) 
has only a pole of the order p + 1 at infinity, i.e. V(x) is a polinomial of degree 
p + 1. Then it remains to consider separately only two different asymptotics: 
small l a and large l a . Of course, since there are many different l a one actualy 
has a vast variety of possibilities: some of l a 's are small and the rest are large. 
The two extreme cases are when all the Z a 's are either small or large, and they 
are refered to as the "character phase" and "Kontsevich phase" respectively. 
The word "phase" is used instead of the more exact "asymptotics" in order 
to | em phasize the relation of these two cases to the strong and weak coupling 
phases in lattice models of Yang-Mills theories. We refer to [35] for some more 
di pcTj ssion of these phases and their properties, here only some basic things will 
be mentioned. 

In the character phase the main observation is that the r.h.s. of the Itzykson- 



Zuber formula (3.6) is essentially the ;hj,racter of the group element g = e of 
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SL(n): 

XR(e L ) = ^^A(x), (3.16) 

where the set {x a } specifies representation R of the SL group 11 . Accordingly 
the integral (3.1) can be represented as a linear combination of characters with 
various R, with coefficients depending on the choice of potential V(x) and N. 
This explains the reason why this limit (when everything is assumed to be 
expandable in positive powers of l a or e La ) is refered to as the character phase. 
One should keep in mind, of course, that the most natural, from this point of 
view, is the situation when integral over x a 's is changed for a discrete sum over 
integer a; a 's - this is one of the possible directions of the search for "quantum 
deformations" of GKM. 

Let us now turn to the limit of large l a . Then the natural expansion would 
be in negative powers of the arguments, but the integral (3.1) does not have 
such expansion, as it is. One should first extract a "quasiclassical factor" and 
it will be the remaining "quantum part" that will possess such an expansion. 
This quantum piece is the most interesting one, and it is for it that the name 
GKM is usualy used. So, in Kontsevich limit, integral (3.1) can be evaluated 
with the help of the steepest descent method, with the classical solution defined 
from V'(Xq) = L. (Note that when doing so we include the logarithmic piece 



11 Correction factor A(l) / A(e l ) can be restored if one considers appropriate generalization 
of Kontsevich integral (which describes not just a puncture but a hole on a surface - in 
terms of the naive string theory), see, for example, ref. [20]. Actually one needs to substitute 
J [dU]exp (tiXULU^ by the loop integral 

[dU(s)] exp (^J ds trX (u(s)a s U^(s) + U(s)LU^(s)) \ = ^p-^ 



When this integral is evaluated for integer x a 's one essentially substitutes every item in the 
product A(Z) = J^ a<t) ('a — lb) by a new infinite product over harmonics, l a — lb — * 
ll k (ia ~h + 2irik) ~ sinh (i^) ~ e l ° - e l " 



in the measure, not in the action.) Let us call solution of this equation Xq = A. 
One could use A from he very beginning instead of L to parametrize Kontsevich 
integral, writing trV'(A)X instead of trLA in the exponent in (3.1). This is a 
natural parameter in Kontsev ich phase, while L = V'(A) plays this role in the 
character phase. The quasiclassical contribution to (3.1), i.e. exponent of the 
classical action divided by determinant of quadratic fluctuations, is equal to: 

■(detA)^ (3.17) 



C V {A\N} = (27T) 1 



exp 


tr| 


AV'(A 


0- 


-V(A))\ 




A 


/detl 


■"( 





and partition function of GKM is by definition 

9 Ml. VI - ^HAjiV} 

ZvAMN}= Cv{Am (3-18) 

This function can be expanded in negative integer powers of A a (i.e. in fractional 
negative powers of original l a ), at least as a formal series. Moreover, there is 
a symmet |gy-| between all the eigenvalues A a , thus Z GKM is in fact a function 
(formal series) of the "time-v? .riables" (the name is historic trace from the theory 
of integrable hierarchies) 

T k = hvA- k (3.19) 

Remarkably, if considered as a function of these T^'s Z becomes independent of 
n! We refer to [4] and [1] for more details (including exact definition of V"(A) 
in (3.17)). 

What remains to be considered here is the eigenvalue representation of Z. If 
(3.12) is divided by the quasiclassical factor Cy{A|iV}, we get: 
■7 ikt T\ 1 det ab y> a+jv (A b ) 

Zv{7V ' T} = (toXF W) — ■ ( } 

Extraction of the quasiclassical factor converts (p(l) into the properly normalized 
expansions in negative integer powers of A: 



e -AV'(A)+V(A) JV"(X\ 

p„(A) = / - V K} <Pa(V'W) = A-Hl + ^(A- 1 )), (3.21) 

V 27T 



It also changes A(7) = A(V'(A)) in the denominator of (3.12) for A (A) ii (J .20). 



Note that, as a corolary of normalization condition (3.21), whenever one puts 
A„ = oo the n x n determinant in (3.20) just turns into 'he (n — 1) x (n — 1) 
determinant of the same form. One can easily understand that this implies the 
n-independence of Z as a function of T- variables (3.19). 

Instead of the simple recurrent relations (3.14) for ip the normalized functions 
ip satisfy: 

<p a (\) = Apa-i(\) = A a -^(X), (3.22) 
where $(A) = ^i(A) and operator . . 

A— 1 9 ^V'"W , A f3 23) 

V"(\) d\ 2{V"{\)f + {6 6) 

now depends on the shape of potential V(x). 

Consider the vector space Tt of all formal Laurent series in some variable A. 
The set of functions 

{A- Q , A"" 1 | a=l,2,..^L 

is one of the many possible basises in this vector space. Let us consider a 
"half-space" H+ consisting of all the series in non-negative powers of A. Then 
{A a_1 | a = 1,2, . . .} is a possible basis in Ti + . Every rotation of the linear 
subspace Ti + within entire Ti can be represented by projection of some basis 
in rotated subspace onto original one. In other words, any semi-infinite set of 
functions {0 a (A)}, such that 

a (A) = A^a + ©(A- 1 )) = A a_1 (l + ]T S ab \~ b ) (3.24) 

f»0 



can be considered as describing some particular rotation of Tl+ in Ti. Of course, 
the same rotation can be represented by different matrices S a b, and in fact ro- 
tations are in one-to-one correspondence with the factor of the set {4> a } modulo 
triangular transformations, which has the natural name of the Universal Grass- 
mannian QTZ. Eq.(3.20) is obviously invariant under such triangular transforma- 
tions, thus Z, as a function of its argument V(x) can be considered as a function 
on QTZ (if the shape of potential is changed, the set {4> a } is also changed). 

Normalization condition (3.24) is, however, invariant not only under tringular 
transformations in {4> a }, but also under the changes A — > A(l + 0(A~ 1 )). Such 
transformations change the point of Grassmannian and they also induce a trian- 
gular linear transformation of time- variables: — ► Tk + lin(Tfc + i, • ■ ■)• 
In other words, Z depends on the choice of variable A on the "spectral curve" 
and on the point of QTZ, i.e. is essentially a function on the tensor product 
QTZ x QTZ of two different Grassmannians. One of them is a space of various 
models (related to the choice of potential in GKM) , another - specifies the basis 
in the space of observables in a given model (related to the choice of time- 
variables). As one expects from the physical arguments and as we just saw 
on a more formal level these two dependencies are in fact interrelated. When 
consideration is restricted to the set of GKM's (from that of all the models 
of string theory) a much more definite statement can be made [36]. Being a 
priori a function of two distinct types of variables: the times {Tf~ = ^trA _fc } 
and potential V(x), the GKM partition function in fact depends only on the 
type of singularities of V(x) (which is a kind of "discrete" information) and on 
peculiar combination of these variables. If V(x) is a polinomial of degree p + 1 
(i.e. has only finite order pole at infinity), then 

Z v {T k } - t p ( jtiW(A)- k / p + P res W^-^dx) , (3.25) 



where W(x) = V'(x), and the shape of the function t p depends only on the 
value of p. (In order to substitute "~" by "=" in (3.25) one should slightly 
redefine the quasiclassical factor and thus Z: one should in fact work with the 
variable L x l p = W(K) x / p instead of A. As often happens, different variables are 
nice for different purposes.) 

If all the arguments Xfc with k > p are put equal to zero, we get a reduced 
r-function 

f p (Ti, . . . ,T p ) = T p (T k )\ Tk= for k>p 

It appears to be a solution to "quasiclassical KP-hierarchy" , which arises from 
pure algebraic construction and can be also identified as partition function of 
topological Landau-Ginsburg model with the superpotential W(x). 



3.6 Kac-Schwarz problem 



The function T p (Tfe) is of course very far from ar Ditrary. First of all it possesses 
peculiar determinant representation of the type (3.20), which in turn implies 
some restrictive bilinear (Hirota) equations and as result r p appears to be a 
KP r-function. Second, r p is further distinguished even among r-functions by 
peculiar features of the functions tp a {X). In particular case of GKM one of 
the ways to represent these properties is to use the recursive relations (3.22), 
<f a = A a ~ 1 Q, suplemented by another obvious property W(A)$(A) = (p p+ i(X) 
(it follows from invariance of the integral ^(A) under the change of integration x- 
variable) . These two relations together give rise to an equation on ipi (A) — ^(A) : 
(A p — W (A))Q(A) — 0, which is just a p-th order differential equation. 

In this way one specifies non-perturbative partition function of some string 
model (in the case of GKM this is in fact a (p, l)-minimal model coupled to 
2c?-gravity) in terms of invariant points of certain operators acting on the Uni- 
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versal Module space Q1Z. This reformulation, though far more abstract than 
the original one, can be very useful for non-trivial generalizations - and be a 
natural step in the search for generic configuration space of the string theory. 
It was first introduced by V.Kac and A.Schwarz [37], and is not yet studied as 
deep as it deserves, even if one deals only with the space of KP r-functions, 
determinant formulas and ordinary Universal Grassmannian. In this (actually 
somewhat narrow) context, the general problem is to describe common invariant 
points in QTZ of two operators, acting on formal Laurent series in A: Va 

Ma G Span{0 b }, (3.26) 



JC4> a eSpan{^} (3.27) 

In the case of GKM A and JC are differential operators of the 1-st and 0-th 
order respectively. Moreover, A is a "gap-one" operator (the gap is equal to 
g if A(f>a — Ylb=i •Aab4'b)- See [19] for some more comments on the gap-one 
case. When both different from unity, the system (3.26)-(3.27) usualy 

describes a multi-parametric set of invariant points, the simplest example being 
associated with (p, g)-minimal models (where p and q are in fact the values of 
gaps for K, and A respectively). In this situation Kontsevich integral desribes 
only duality transformation of (p, g)-model into the (q,p)-one [38] (note that 
these models do not coincide after they are coupled to 2d gravity). 



3.7 Ward identities for GKM 

Advantage of the Kac-Schwarz reformulation of GKM is that it is very easy to 
deform, since there are no real restriction imposed on the choice of operators 
A and K,. However, instead, this formulation does not introduce any reach 
structure and does not immediately provide any valuable information about the 



form and properties of solutions to (3.26)-(3.27), i.e. does not explicitly reveal 
any nice properties which could be common for all the string models. Thus it 
could serve as a starting, but not the final point of analysis of non-perturbative 
partition functions. Alternative approach, making use of the Ward identites, 
provides a better description of GKM, but can appear too restrictive to allow 
for any interesting deformations. From the point of view of Grassmannian, the 
Ward identities specify some subset of points in Q1Z, which are invariant under 
certain subalgebras oilAGL(oo) (the symmetry group of entire Grassmannian). 
The corresponding "homogeneus spaces" are normally not discrete and contain 
a vast variety of points. By themselves the Ward identites are not enough to 
specify particular points in Q1Z uniquely. They should still be supplemented 
by some extra conditions, like "reduction constraints" (in fact one should keep 
one of the Kac-Schwarz constraints, (3.27), and only another one can be usualy 
substituted by the symmetry- like relation, coming from the Ward identities). 



3.7.1 Gross-Newmann equation 

We refer to [1] for a very detailed review ^FJhe Ward identities in GKM. They 
are all corollaries of the simple Gross-Newjnann (GN) equation [39], imposed 
on Kontsevich integral (3.1) as result of tlsj invariance under arbitrary change 
of the integration variable X : 

J d Xe^ v ^ +m ^ x+t - [LX (-V'(X) + NX- 1 +L)=0. 

or 

(V {d/dLtr) -L-N {d/dL tr )~ l \ T v = 0. (3.28) 

Being just equation of motion for !Fv{L} eq.(3.28) provides complete informa- 
tion about this function. However, this statement needs to be formulated more 
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carefully. One of the reasons is that (3.28) does not account explicitly for a very 
import ant property of that it actualy depends only on eigenvalues of L. 

This information should still be taken into account explicitly. If this is kept in 
mind, it becomes a tedious but straightforward work to substitute Ty — CyZy 
and express L-derivatives through T^-derivatives in order to derive Virasoro and 
ly-constraints in conventional form [1]. We shall briefly sketch some pieces of 
this derivation and related problems in the remaining part of this section. 



3.7.2 VT-operators in Kontsevich models 

First of all, the Gross-Newmann equation (3.28) for Kontsevich models can be 
easily expressed in terms of tli^-sp called VK-operators. Namely, we shall prove 
the following identity [40] : 

/ a x m+1 

I ° \ VST. \ - A-i-ln/(™+ 1 ) 



\dA t 



Z{T k } = (±) m+1 ^A-'- 1 H> ; ( ";t 1) ( r )2{T fc }, (3.29) 



z>o 



valid for any function Z which depends on T k = =p^trA j% j k > 1 and 
T = itrlogA. 

The H^joperators are defined [40] by the following construction. Consider the 
action of Tr^^-L" on e T ' cU ^ L ' , = * fcTri . It gives some linear combination 
of terms like 



trZ 0l ...tr£° ! e tr£ C£f = n n e~ tlU ^ (3.30) 
dt ai ...dt ai 

i.e. we obtain a combination of differential operators with i-derivatives, to be 
denoted W(t): ^ 

Wti\t)e tru ^ = Tr^—L n e w ( L \ m,n>0. (3.31) 
oL tr 



For example, 



see the next subsection. In other cases calculations with the use of identity 



WW = J±- n > o- 

dtn' 



(3.29), accounting for the quasiclassical factor Cv{L} and_t 



£ 
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ic difference be- 



, n>-l; 



fe=0 



^n 3) = £ ktkHl 



k,l = l 



dtk+n dtkdt n -k 

K—{) 

d °° 

— hVfctfc V - 

dt k+ i +n f-J d< Q di h 



<9 2 



■£>* £ 



, dt a dt b+k _ 1 ' ^ dt a dt b dt c 

k—1 a-\-b—n-\-l a+o+c— n 



E 



c> 3 



(n + l)(ra + 2) 9 
2 



(3.32) 

Note, that while Wn and Wn are just the ordinary (/7(1)-Kac Moody and 
Virasoro operators respectively, the higher W-^" 1 ) -operators do not coincide with 
the generators of the W- algebras: already 



w?^w^= E^7^ + 2 E fc ^ E " 



, , . dtk+i+n ~t . . dt a dt b 

k,l=l k—1 a-\-b—k-\-n 



4 E 



9 3 



(3.33) 



3 , It"' dt a dt b dt c 

a-\-o-\-c—n 



^-operators (in variance with ordinary PF-operators) satisfy recurrent relation: 

m+n— 1 

, . ,-, (3.34) 



k=l k=0 k 



Actually not too much is already known about the W operators and the struc- 
ture of W-algebras (in particular it remains unclear whether the negative har- 
monics Wn m+1 ^ with n < —m can be introduced in any reasonable way), see 
[40] for some preliminary results. 

Now we can come back to the identity (3.29). Its most straightforward 
application is to the Gaussian Kontsevich model with potential V{x) = \, 



tween L = V'(A) and A become somewhat more involved, though still seem 
sufficiently straightforward. Also for particular potentials V(X) partition func- 



tion Zy{T} is a 



ally independent of certain (combinations of) time-variables 



(for example, if V(X) 



x p+1 
P+i > 



it is independent of all the T p k, k G Z + ), and 



this is important for appear enge of the constraints in the standard form, i.e. for 



sttaJr 



certain reduction of H^-consttaihts to the ordinary VK-constraints. This relation 
between W- and W-operators deserves further investigation. 

The proof of eq.(3.29) is provided by the following trick. Let us make a sort 
of Fourier transformation: 



Z{T} = JdH g{H}e^=o TkTtH \ 



(3.35) 



where integral is over N x N Hermitean matrix H . 12 Then it is clear that once 
the identity (3.29) is established for Z{T} substituted by e TrU[H \ U(H) = 
SfeLo ^fcTrif fe , with any matrix H, it is valid for any function Z{T}. The 
advantage of such substitution is that we can now make use 



of 



the definition 



(3.31) of the W operators in order to rewrite (3.29) in a very explicit form: 

m+l 00 

s ™ ( tf) =(±r+ l£ A -*-l^_ m 



d 



dA 



tr 



00 

-Ea-'- 1 

2>0 



W^ 1] {T)e TrU[ - H) 



(±) 



OO 

E-\ 

l>0 



1-1 



Tr 



d 



dH tr 



H l e TrU(H) 



(3.36) 



(±) 



m+l 



Tr 



d 



tr 



1 



A® I — I <g) H 



12 Here it is for the first time that we encounter an important idea: matrix models - the 
ordinary 1-matrix model (2.1) in this case - can be considered as defining integral trans- 
formations. This view on matrix models can to large extent define their role in the future 
developement of string theory. 
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Now expression for T's in terms of A should be used. Then 



e TrU(H) = Def ±l( A 8/ _/ 8ff ) 



(3.37) 



and substituting this into (3.36) we see that (3.29) is equivalent to 



dA tr 



m+l 



(±) 



m+l 



7-Tr 



( 9 



1 



Det ±l (A®I-I®H) 



A® 7- 7® 77 




(3.38) 



Here "Tr" stands for the trace in the 7f-space only, while Det = Det ® det - for 
determinant in both H and A spaces. After one A-derivative is taken explicitly, 
we get: 



(7 ® Tr) 



( d 



I-Ii 



\dA tr 

Det ±:L (A® I- J® 77 
A®7-7®77 



±- 



d 



dH tr 
= 0. 



(3.39) 



This is already a matrix identity, valid for any A and H of the sizes n x n and 
N x N respectively. For example, if m = (M^^-casc), |t~i| obviously satisfied. 
If both n = N = 1, it is also trivially true, though for different reasons for 
different choice of signs: for the upper signs, the ratio at the l.h.s. is just unity 
and all derivatives vansih; for the lower signs we have: 



d_ 

ox 



d_ 

~~dh 



( 



y (-VY-- 



\ a,b>0 



dh 



J 



and this obviously vanishes since (^ + ^j)/(A — h) = for any f(x). 

If rn > and A, H are indeed matrices, direct evaluation becomes much more 
sophisticated. We present the first two nontrivial examples: m — 1 and to = 2. 



The following relations will be usefull. Let Q = j^gj^j^jj- Then 
Det^Q^—Det^Q = ± [(7® Tr)Q] ; 

OJXtr 

Det^Q-^—Det^Q = =p [(tr ® I)Q] ; 

Otitr 

° I)Q = - [(tr ® I)Q] Q; 



dA tr 

d 



(3.41) 



,)77 6r 



= [(7 ® Tr)Q] . 



This is already enough for the proof in the case of to = 1. Indeed: 
Det^Q ( ® I T I ® J^A QDet Tl Q = 

V dA tr dH tr J 

= {- [(tr ® I)Q] Q ± [(I ® Tr)Q] Q}t 



(3.42) 



T{[(7 ® Tr)Q] Q T [(tr ® I)Q] Q} = 0. 
The first two terms at the r.h.s. come from A-, while the last two - from 
77-derivatives. 

In the case of to = 2 one should take derivatives once again. This is a little 
more tricky, and the same compact notation are not sufficient. In addition to 
(3.41) we now need: 



( 9 



\dAt 



I I [(tr ® I)Q] Q = - [(tr ® I)Q] 2 Q-B. (3.43) 



Here 

[(tr ® I)Qf - [(tr ® I) [(tr ® I)Q] Q] , (3.44) 

while in order to write B explicitly we need to restore matrix indices (Greek 
for the A-sector and Latin - for the H one). The (ai, 7fc)-component of (3.43) 
looks like: 

{dxy/" 1 ) Q ™s Q & = - Q M l ?Q% - q%q«sQ% (3-45) 



and appearence of the second term at the r.h.s. implies, that B^L 
Q%Q l LQ%- Further, 

/ ^ i) [{I ® Tr)Q] Q = 



\dA t 

[(I <g> Tr) [(tr ® I)Q] Q]Q - [{I <g> Tr) [(/ ® Tr)5tt] Q; 



[(tr®7)Q]Q = (3.46) 



377 tr , 

+ [(tr <8 I) [(I <8 Tr)Q] Q] Q + [(/ <8 Tr) [(tr ® 7)Q] Q] Q; 

5 \ 

® iF" ® Tr )Q] Q = + P ® Tr) [(/ ® Tr)g] Q Q + 6. 
an tr J 

It is important that ,B that appears in the last relation in the form of = 
QagQfpQp-y i s exactly the same £? as in eq.(3.43). 
Now we can prove (3.39) for m = 2: 

Det±lQ {i^;) 2 ® 1 - 1 ® {ok)) QDetTlQ = 

= {± [(/ ® Tr)Q] (- [(tr <g> I)Q] Q ± [(/ ® Tr)Q] Q) - 

- (- [(tr ® 7) [(tr ® I)Q] Q]Q-B)± 

± (- [(I ® Tr) [(tr <8> 7)Q] Q] Q - [(tr® 7) [(7® Tr)Q]Q]Q)}- ( -' 1 ' 4 ' ' 

- {T [(tr ® 7)Q] ([(/ (8 Tr)Q] Q T [(tr <8 7)Q] Q) + 

+ ([(/ ® Tr) [(/ (8 Tr)Q] Q] Q + 6) T 

T ([(tr <g> /) [(/ <g> Tr)Q] Q]Q + [(I ® Tr) [(tr <g> 7)Q] Q] Q)} 
where the terms 1,2,3,4,5,6 in the first braces cancel the terms 1,3,2,4,6,5 in the 
second braces and identity (3.44) and its counterpart with (tr ® 7) — > (7 <g) Tr) 
is used. 

Explicit proof of eq.(3.39) for generic m is unknown. 
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3.7.3 Discrete Virasoro constraints for the Gaussian Kontsevich 
model 



As a simplest illustration we derive now the constraints for the Gaussian Kont- 



sevich model [41] with potential V(X) — IX 



Z^{N,T} = 



dX(detX) N e 



{detL) N 

In this case L = V'(A) = A, and the time- variables are just 



(3.48) 



k 



k 



(3.49) 



The model is non-trivial becuase of the presence of "zero-time" variable N [42]. 
The Gross-Newmann equation (3.28) looks like 



( 9 



n+l 







- N 



d 



□ 



(detL) N \dL tr J \dL tr " \dL tr 
■{detL) N e +tI ^Zx2{N,T} = 



-L • 



(3.50) 



In order to get rid of the integral operator (^) _1 one should take here n > 
rather than n > — 1. In fact all the equations with n > follow from the one 
with n — 0, and we restrict our consideration to the la st o ne. For n — we 
obtain from (3.50): 



r d n 



\\dLt 



2N -L 



f d N 
\dL tr L 



□ 

□ 

L ) \Z = (3.51) 



One can now use eq.(3.29) to obtain: 



E r,m+2 E ( tr r,fc 



1 \ d 



E 



L m+2 i ^ V L k J dT k+m ' ^ dT k dT m+k 

d d / 1 \ \ 

— 2iV-— + N 2 S m . Q - JV tr- S m ,-! Z = 

dT m+2 dT m \ LJ 

= E 7 ^e^X m (T + r)e-^ Z = . 



(3.53) 



m— — 1 

K2) 



Here L m (t) = Wm (t) are just the generators (2.9) of "discrete" Virasoro algebra 
(2.8): 




- Nt0 . (3.54) 



and at the r.h.s. of (3.53) r k = — |#fe,2- 13 

Thus we found that the Ward identities for the Gaussian Kontsevich model 
(3.48) coincide with those for the ordinary 1-matrix model (2.1), moreover the 
size of the matrix N in the latter model is associated with the "zero-time" in the 
former one. This result [41] of course implies, that the two models are ide itid al: 



e- NT °Zx2{N,T u T 2 , . . .} ~ Z N {T , T l7 T 2 , . . .}. 



See [42] and [1] for explicit proof of this identity. 



(3.55) 



This small correction is manifestation of a very general phenomenon which was already 



mentioned in s.3.5 above: from the point of vie^i 



natural to consider Zy not as a function of T-v ariab les. but of some more complicated com 



Hi symmetries (Ward identities) it is more 



bination Tj, depending on the shape of potential V. If V is a polinomial of degree p+ 1, 

T k = £tr(V (\))~ > p , while r k = k ^_ k ^ Res (V'(fj,)) " dfJ,. For monomial potentials these 
expressions become very simple: T k = T k and r k = — jq^y<5fc, p +i. See [36] for more details. In 
most places in these notes we prefer to use invariant potential-independent times T k , instead 
of T k , but then Ward identites acquire some extra terms with r k . 
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x 3 

3.7.4 Continuous Virasoro constraints for the V = 4p Kontsevich 
model 

This example is a little more complicated than that in the previous subsection, 
and we do not present calculations in full details (see [43] and [4]). Our goal 
is to demonstrate that the constraints which arise in this model, though still 
form (Borel subalgebra of) some Virasoro algebra, are different from (2.8). 
From the point of view of the CFT-formulation the relevant model is that of 
the twisted (in this particular case - antipcriodic) free fields. These so called 
" continuous Virasoro constraints" give the simplest illustration of the difference 
between discrete and continuous matrix models: this is essentially the difference 
between " homogeneous" (Kac-Frenkel) and " principal" (soliton vertex operator) 
representation of the level k = 1 Kac-Moody algebra. From the point of view of 
integrable hierarclrrel this is the difference between Toda-ch|iIn|-like and KP-like 
hierarchies. 

Another (historicaly first) aspect of the same relation also deserves mention- 
ing, since it also illustrates the interrelation between different models. The 
discrete 1-matrix model arises naturally in description of quantum 2d gravity 
as sum dvnr 2-geometries in the formalism of random equilateral triangulations. 
The model, however, decribes only lattice approximation to 2d gravity and 
(double-scaling) coninuum limit should be taken in order to obtain the real 
(continuous) theory of 2rf |-g¥p vity. This limit was originally formulated [44] in 
terms of tl e contraint algebra (equations of motion or " loop" or " Schwinger- 



Dyson" equations - terminology is taste-dependent), leaving open the problem 
of what is the form of partition fnncti fm-g 00 "!"^! of continuous theory. Since 
the relevant algebra appeared to be just the set of Ward identities for Kontse- 



vich model (with V(X) 



□ 



this pjrojes that the latter one is exactly the 
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continuous theory of pure 2c? gravity. At the same time, Kontsevich model itself 
can be naturally introduced as a theory of topological gravity (in fact this is 
how the model was originally discovered in [3]). From this point of hfWw the 
constraint algebra, to be discussed below in this subsection, plays thfe-etentral 
role in the proof of equivalence between pure 2d quantum gravity and pare 
topological gravity (in both cases "pure" means that "matter" fields arc'-not 
included) . 

After these introductory remarks we proceed to calculations. Actually they 
just repeat those from the previous subsection for the Gaussian model, but 
formulas get somewhat more complicated. This time we do not include zero- 
time N and use eq.(3.28) with V(X) — 4p. Also, this time it is much more 
tricky (though possible) to work in matrix notations (because fractional powers 
of L will be involved) and we rewrite everything in terms of the eigenvalues of 
L, 



Substitute 



i-r 2 \ 3/2 



( d 2 \ d 2 , ^ 1 ( 9 d 



(3.56) 



and introduce a special notation for 

V -r- 1 9 c 9 + /r 1 1 V 



VX a "^dXa""? dX a ' v " a 4A a 2^ VKiVh + VK)' 



(3.57) 



Then (3.28) turns into 



' " ^'ET-^:f^:-^:))^m = o. <*»> 



\VXaJ ~ X a — Xb \VX a VX b J 



Now we need explicit expression for T: 

T k = jL~ k , (3.59) 
k 

and - as we already know from the previous subsection - we also need 



n = -3^,3. (3.60) 

2x±{T} is in fact independent of all the time- variables with even numbers 
3 

(subscripts), see [4], [1] for the explanation. Therefore we can take only k = 
21 + 1 in (3.59), 



T21+1 = 



21 



(3.61) 



r 2 i+i = ~3^.l 



and 



d 2 



9 ? lr\ V- ^+i dZ l^.-i-* dZ 



9A2 



dX a 



d T 2l+1 dZ _1 
dX a dT 2l+l 2 

1=0 a M+L a=0 

dZ 
dT 2l+1 dT 2 

m+1 ^ - - * 



dT 2l+1 



l—m—S 



/=o 



dT 2l+1 ■ 
(3.62) 

These expressions should be now substituted into (3.58) and we obtain: 

dZ 



m— 3 



l.m=0 



□ 



1=0 
00 



-1-* 

2 \ : 



a=0 



b=£a 



1 1 



(3.63) 



dZ 



with 



coefficient, reordering the items: 



C 2n = E ( 1 + o ) ( T2i + 1 + r2i + 1 ) 



d 



d 2 



i l 2r 



a i+ ^ n i dT 2l+l dT 2m+l i6 r 

Z,m>0 



1 oo 



d 1 



a 2 i l. 



odd fc=l odd fc— 1 

(3.64) 

Factor | in front of the first term at the r.h.s. in (3.64) is important for C 2n to 
satisfy the properly normalized Virasoro algebra: 14 



[^2nj^2m] — (fl — m)£ 2 n+ 2 m- 



(3.65) 



Coefficient j in front of the second term can be eliminated by rescaling of time- 
variables: T — > IjT, then the last term turns into j*T^5 n -i. 

We shall not actually discuss evaluation of the coefficient in front of Z (with 
no derivatives), which is denoted by [. . .] in (3.63) (see [43] and [4]). In IfaJtr— | 



almost all the terms in original complicated expression cancel, giving finally 



T 2 



(3.66) 



L "' J 16A2 4A a ' 

and this is represented by the terms with 8 n fi and S n ^-i in expressions (3.64) 
for the Virasoro generators C 2n . 

The term with the double T-derivative in (3.63) is already of the necessary 
form. Of intermidiate complexity is evaluation of the coefficient in front of 



dZ 



>)T 2 ; 



in (3.63), which we shall briefly describe now. First of all, rewrite this 



14 Therefore it could be reasonable to use a different notation: C n instead of C2n ■ We prefer 
£271 > because it emphasises the property of the model to be 2-reduction of KP hierarchy (to 
KdV). 



b^a 



4 Aa 



\E 



-1-2 



2^ VXb 



(3.67) 



- A 



-1-1 



The first two terms together are equal to the sui ji oj ver all b (including b = a): 



£ 



2 A a — Ah 

b 



-l-S. -1-2 
A a 2 - A ' 2 



2\' a 



1 /4- 3 

I \ /+ ! \ 2 \ 5 

^ A a f H a b 

1+2 



X 



1 \ L ^ 2 v ' : 

1 v A a — A Q 

2 V ^a-Xb 
1 



\ l ' 2 \ 1 ' s 

Aa A, 



A,, 



A 



(3.68) 



Similarly, the next two terms can be rewritten as 

A a 

6 



9 ZLJ/r. 



a q + vAb 2A Q 



1+2 £ ' 



A„ — A(, 



□ 



2A a ~ 



i+2 £ 



\ 2 \ l+ 2 \ ( 
Aa A,, — A. 



(3.69) 



Aa — Ah 



l+z 

\ * 2 



The sum of these two expressions is equal to 



2X l a 



1 \ -» 

1+2 Z^ 



Ai+ 2 - A<+ 2 1 



A a — Ah 



A 



(3.70) 



b « 6 

Note that powers ^ + 2 are already integer and the remaining ratio can be 
represented as a sum of / + 2 terms. Adding also the last term from the l.h.s. 
of (3.67), we finally obtain: 



Ai+ 



!_ if _Lv— L_ 

b 



9 \n 

1 n=-l Aa 



b A 



\ £ Ti2(2a-2n+l)(T + r) 2i _ 2 „ +1 

z , A a 



(3.71) 
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in accordance with (3.63) and (3.64). 

Calculations can be repeated for every particular monomial potential V{x) — 
^pj-, but they become far more tedious and no general derivation of W^- 
constraints [44] is yet found on these lines. See [45] for detailed examination of 
the -constraints in the -^--Kontsevich model. 



4 KP/Toda r- function in terms of free fermions 



There are several different definitions of r-functions, but all of them are par- 
ticular realizations of the following idea: T-function is a generating functional 
of all the matrix elements of some group element in particular representation. 
Since methods of geometrical quantization allow to express all the group the- 
oretical objects in terms of quantum theory of free fields, generic r-functions 
can be also considered as non-perturbative partition functions of such models. 
The basic property of T-function, which can be pracJteJlly derived in such a 
general context, is that it always satisfy certain bilinear equations, of which 
Hirota equation for conventional KP T-function is the simplest example. 

with the free particles of a peculiar 
~[4o]. Existence of fermionization is 



KP/Toda t- functions are associated 
type: free fermions in 1 + 1 dimensions" 
a very rare property of free field theory (in varience with bosonization which is 
always available). If existing it leads t^~^lramatic simplification of the for- 
malism and to especially simple determinant formlulas (instead of sophisti- 
cated and often somewhat abstract objects like chir ^rT^ eterminants det d in 
generic case). In the case of Kac-Moody algebras the corresponding T-function 
is nothing but non-perturbative partition function of the corresponding Wess- 



Zumino-Novikov-Witten model. Among simply-laced algebras only GL{N) k=1 
is straightforwardly fermionized, and the formalism is much simpler in this 
case than for generic Wess-Zumino-Witten model with arbitrary level k. For 
N = 1 we obtain KP/Toda T-functions, while N ^ 1 are related to the "N- 
component KP/Toda systems". Level-one Kac-Moody algebras SL(N) k=1 are 
distinguished because their universal envelopping are essentially the same as 
those of their Cartan subalgebras. This allows to define generation functions 
with the help of sets of mutualy commuting generators and makes evolution, 
described by commuting Hamiltonian flows, complete (acting transitively on the 
orbits of the group) . This is why such systems are distinguished from the point 
of view of Hamiltonian integrability - and why they are the usual personages in 
the theory of integrable hierarchies. In general case (k ^ 1) one naturally deals 
with the set of flows that form closed but non-Abelian algebra. In the language 
of matrix models restriction to k = 1 and free fermions is essentially equiva- 
lent to restriction to eigenvalue models. Serious consideration of non-eigenvalue 
models, aimed at revealing their integrable (solvable) structure will certainly 
involve the theory of generic T-functions. 

□ 

4.1 Explicit definition 

Let us introduce two fields (a spin-1/2 b, c-system) ip(z) and tp(z) satisfying 
canonical commutation relation: 

ip(z)]+ = S(z - z)dz 1 / 2 dz 1 ' 2 . (4.1) 

Then 

t{A} ~ (0 | exp ( I I A(z, ~z)i/j(z)4>(~z)) | 0). (4.2) 

\Jdz Jdz / 



Now it is usual to expand in Laurent series around z = 0: 

n&Z nEZ 
\^prm '0n]+ — $m,n 5 

V>m I 0) = for m < 0; ^ m | 0) = for m > 0; 
A(z,z)= 51 «~ m_1 5 n >4mnd« 1/a d2 :l/a ; 



(4.3) 



m,n£Z 



so that 



A(z, Z)i>(z)i>(z) = ^ A mni>mi>n- (4.4) 
diJ(fa m.neZ 

In fact this expansion could be around any pair of points zq, z^ and on a 
2-surface of any topology: topological effects can be easily included as specific 
shifts of the functional A(z, z) - by combinations of the " hadle-gluing operators" . 
Analogous shifts can imitate the change of basic functions z n for z n+a and more 
complicated expressions (holomorphic 1/2-differentials with various boundary 
conditions on surfaces of various topologies). 

One can now wonder, whether local functionals A(z,z) = U(z)S(z — 
z)dz 1 / 2 dz 1 / 2 play any special role. The corresponding contribution to the 
Hamiltonian looks like 



Hcartan = f U {z)^(z)^( Z ) = f U(z)J{ Z ), (4.5) 

where 

J(z) = ip(z)ip(z) = ^2 J n z~ n ~ 1 dz (4 g) 
is the U(l)k=i Kac-Moody current; 

Jn = i ) m'ipm+n\ \Jmi Jn] = m^m+nfi- (4-7) 

mez 



□ 



If scalar function (potential) U(z) is expanded as U(z) = J^kez ^kZ k , then 

□ 

Hcartan = tkJk- (4.8) 

n£Z 

□ 

This contribution to thC3hole Hamiltonian can be considered distinguished for 
the following reason. Let us return to original expression (4.2) and try to con- 
sider it as a generating functional for all the correlation functions of ip and tp. 
Naively, variation w.r.to A[z, z) should produce bilinear combination ip(z)if)(z) 
and this would solve the problem. However, things are not just so trivial, be- 
cause operators involved do not commute (and in particular, the exponential 
operator in (4.2) should still be defined less symbolically, see next subsection). 
Things would be much simpler, if we can consider commuting set of operators: 
this is where abelian U(l) k=1 subgroup of the entire GL{oo)k=i (and even its 
purely commuting Borel subalgebra) enters the game. Remarkably, it is suffi- 
cient to deal with this abelian subgroup in order to reproduce all the correlation 
functions. 15 The crucial point is the identity for free fermions (generalizable to 
any b, c-systems): 



: V(AW(A) : = : 




5 We once again emphasize that this trick is specific for the free fermions and for the level 
k = 1 Kac-Moody algebras, which can be expressed entirely in terms of free fields, associated 
with Cartan generators (modulo some unpleasant details, related to "cocycle factors" in the 
Frenkel-Kac representations [12], which are in fact reminiscents of free fields associated with 
the non-Cartan generators (parafermions) [47], - but can, however, be put under the carpet 
or/and taken into account "by hands" as "unpleasant but non-essential(?) sophistications). 



which is widely known in the form of bosonization formulas: 16 if J(z) = d(f>(z 
i>(\) ~ : e^ X) : ( : ^(oo)^(A) : = : e { ^-«"» : \ ; 
^(A) ~ : e"* (A) : ( : ^(A)^(oo) : = : e W°°)-4>W) : ) . 



(4.11) 



This identity implies that one can generate any bilinear combinations of tp- 
operators by variation of potential U(z) only, moreover this variation should be 
of specific form: 




keZ / Jz k£Z Jz 



k k \ z k z k 
k£Z v 



(4.12) 



i.e. 



**> = Kin i) ^ 

Note that this is not an infinitesimal variation and that it has exactly the form, 
consistent with Miwa parametrization. 

Since any bilinear combination can be generated in this way from U(z), it 
is clear that the entire Hamiltonian ^ ^4mn'0m'0n can be also considered as 
resulting from some transformation of V (i.e. of " time- variables" tfc). In other 
words, 

T{A} = A [t]T{A = U}. (4.14) 



le Formulas in brackets are indeed correct, they are preceded by the usual symbolic relations. 
Using these formulas we get: 

: VWCA) : = : : : jl ' : ( 4 - 10 ) 

This identity can be of course obtained within fermionic theory, one should only take into 
account that i/>-operators are nilpotent, so that exponent of a single i/>-operator would be just 
a sum of two terms (polinomial) and carefully follow the normal oredering prescription. 



□ 



These operators O a are naturally interpreted as elements of the group GL(oo), 
acting on the Universal Grassmannian GR [48]- [50], parametrized by the ma- 
trices A mn modulo changes of coordinates z — > f(z). This representation for 
r{yl} is, however, not very convenient, and usually one considers infinitesimal 
version of the transformation, which just shifts A 



r{t\A + SA} = d SA [t}r{t | A}, (4.15) 



note that this transformation clearly distinguishes between the dependencies of 
ront and on all other components of A. The possibility of such representation 
with the privileged role of Cartan generators is the origin of all simplifications, 
arising in the case of free-fermion r-functions. Relation (4.15) is the basis of 
the orbit interpretation of r-functions [49] . 



4.2 Basic determinant formula for the free-fermion corre- 
lator 

Let us consider the following matrix element: 

r N {t,i\G} = (N\e H Ge R \ N) (4.16) 



where 



G = CXp ^2 A mn^ m 1pn 5 
\m,neZ J 

h = *^2tkJk, h — tkJ-k 

fe>0 fe>0 

J(z) = ip(z)-ip(z) = ^2 J n z~ n ~ 1 dz; J n = ^ ^k^k+n] 
neZ k 



ipm I N) = 0, m< N 
ip m | N) = 0, m > N 
J m \ N) = 0, m > 



(7V|V m = 0, m>7V; 
(iV | {p m = 0, m < N; 
(N\ J m = 0, m < 0. 



(4.17) 



The "iV-th vacuum" | N) is defined as the Dirac sea, filled up to the level N: 



N-l 



| N) = Y[ & | oo) = JJ i>i | -co) ; 

z— AT i— — oo 

oo JV-1 

(Ar |= (oo | JJ = (-oo | |"{ 4>i, 



(4.18) 



»=JV 



where the "empty" (bare) and "completely filled" vacua are defined so that: 



tp m | -oo) = 0, (-oo | ip m = 0, 

1pm | OO) = 0, (OO | 1pm = 



(4.19) 



for any m e Z. For the only reason that operators J, iJ, H and G are defined 
so that they always have tp at the very right and tp at the very left, we get also: 

J m | -oo) = 0, (-oo | J m = 0, 

G ±x | -oo) =| -oo); (-oo | G ±x = (-00 |; (4.20) 

e ±H I -00) =| -00); (-00 I e ±H = (-00 I . 
Now we can use all these formulas to rewrite our original correlator (4.16) as: 

(N I e H G e n \ N) = 

= (-00 I ( n i>i J e ff G e B ( JJ ) I -00) = 

\i— — 00 / \«— — 00 / 

= (-00 I e- H ( Y[ ^ e H Ge*(\[ A e~ S | -00) = 

\i— — 00 / \i— — 00 / 

AT-1 AT-1 

= (-00 1 n **[*] n 1 -°°) = 



(4.21) 



= Dct_ 00<iJ<A r(-oo I *i[t]*f [i] I -00) = 

= Detij^Hi+Nj+N- 
The last two steps here were introduction of "GL(oo)-rotated" fermions, 

= e-'tyie*; ^[i] = e % e - S ; *f [i] = G^G" 1 , (4.22) 

and application of the Wick theorem to express multifermion correlation func- 
tion through pair correlators 



Wy(t,t) = (-00 I *i[t]*f [t] I -OO) = 



(4.23) 

= (-00 I *<[t] G *,•[«] I -00), 



(once again the fact that G 1 | —00) =| —00) was used). The only non-trivial 
dynamical information entered through applicability of the Wick theorem, and 



for that it was crucial that all the operators e H , e H , G are quadratic expo- 
nents, i.e. can only modify the shape of the propagator, but do not destroy the 
quadratic form of the action (fields remain free). This is exactly equivalent to 
the statement that "Heisenberg" operators are just "rotations" of ip, i.e. 
that transformations (4.22) are linear. 

We shall now describe these transformations in a little more explicit form. 
Namely, their entire time-dependence can be encoded in terms of the ordinary 
Shur polinomials P n (t). These are defined to have a very simple generating 
function (which we already encountered many times in the theory of matrix 
models) : 

^P n (t)z n =exp(jrt k z k ) (4.24) 

n>0 \fe=l / 

i 2 

(i.e. Po — 1, P\ — t\, Pi = -k + t2 etc.), and satisfy the relation 

dP n „ 



dt k 

Since 



-k ■ 



ex P (f>*W(E^^ 

\fc=l / fe>0 \n k >0 K ' t 



(4.25) 



(4.26) 



Shur polinomials can be also represented as 



Ekni, —n 
k>Q K 

Now, since 



□ 



e- B Ae B = A + [A, B] + ± [[A, B},B} + 1 [[[A, B},B],B} + ... (4.28) 

and 

[4>i, Jk] = -tpi+k, [[i>i,Jk 1 ],Jk 2 ]=i>i+ki+k2i ■•■■> ( 4 - 29 ) 



we have for every fixed k: 



t 

n k >0 



(4.30) 



It remains to note that all the harmonics of J in H = X)fc>o t&Jk commute with 
each other, to obtain: 

*,(*) = e- H ^e H = ( ft e-^* j & [\{e^A = 
\fe>o / \fe>0 / 



n>0 



+71 



{n fc } \fe>0 

V > knu—n 



(4.27) 



/ 



Similarly, relation [J^,^] = implies, that 



ri>0 



and finally 



W iJ = E< -00 I ^ G ^™ I -°°) Pl-i(t)P m -3(t) 



= 5^T lrn fl_ i (t)P m _ i (t), 



which implies also that 



aw. 



1J -W 



— Ti-ij+k- 



Matrix 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



Ti m = (-oo I ipi G ip m I -co) 



(4.35) 



u 



D 

is the one which defines fermion rotations under the action of GL(oo)-group 



element G: 

Gi/imG^ 1 = ^IplTim, 
lez 

G^lG = Tlrn^m, OT G^G" 1 = ^ (T' 1 )^ 

mez mez 
If G = 1, T lm = 8 lm . If all t k = i k = 0, Hij = T iy 



(4.36) 

□ 



□ 



4.3 KP hierarchy and other reductions 

In the previous subsection a formula 

r N {t, t | G} = Det itj<0 H i+ N,j+N ( 4 - 37 ) 

was derived for the basic correlator, which defines " Toda-lattice r- function". 
For obvious reasons the variables t are often refered to as "negative-times" . 
T-function can be normalized by division over the same quantity with all the 
time- variables vanishing, but this is not always convenient. Eq.(4.37) has gener- 
alizations - when similar matrix elements in a multifermion system is considered 
- this leads to "multicomponcnt Toda" (or AKNS) r-functions. Generaliza- 
tions to arbitrary conformal models should be considered as well. It has also 
particular "reductions", of which the most important are: KP (Kadomtsev- 
Petviashvili) , forced (semi- infinite) and Toda-chain r- functions. This is the 
subject to be discussed in this subsection. 

Idea of linear reduction is that the form of operator G, or, what is the same, 
of the matrix T\ m in eq.(4.33), can be adjusted in such a way, that tat{£, i \ G} 
becomes independent of some variables, i.e. equation(s) 



(E a ^r k + E «^ + E 0kD N (k) + 7 j m{t, 1 1 g} = o 



(4.38) 



can be solved as equations for G for all the values of t, t and N at once. (In 
(4.38) Djsi{k)fM = Jn+u — In-) In this case the system of integrable equations 
(hierarchy), arising from Hirota equation for r, gets reduced and one usually 



speaks about "reduced hierf ,rc ly" . Usually equation (4.38) is imposed directly 



on matrix Tiij, of course than (4.38) is just a corollary. 



We shall rkfjdr to the situation when (4.38) is fulfilled for any t,t,N as to 
"strong reduction". It is often reasonable to consider also "weak reductions", 
when (4.38) is satisfied on particular infinite-dimensional hyperplanes in the 
space of time- variables. Weak reduction is usually a property of entire r- 
function as well, but not expressible in the from of a local linear equation, 
satisfied identicall for all values of t, i, N. Now we proceed to concrete exam- 

□ 

Toda-chain hierardiy. This is a strong reduction. The corresponding con- 
straint (4.38) is just 

dt k ~ di k ' l4 -^ j 
or, because of (4.34), Tti+k.j = "Hi.j+k- It has an obvious solution: 

H id =H i+j , (4.40) 

i.e. TLij is expressed in terms of a one-index quantit |y~| ii . It is, however, not 
enough to say, what are restrictions on Hij - they should be fulfiled for all t 
and t at once, i.e. should be resolvable as equations for T; m . In the case under 
consideration this is simple: Tj m should be such that 

Ti m = f l+m . (4.41) 



□ 



Indeed, then 

l,m l,m 

/ „ \ (4-42) 

n>n \fc=n / I — I 



and 



Hi = J2 T n+i J2 Pk i^n-k (t) ■ (4.43) 



n>0 \k=0 



Volterra hierarchy. Toda-chain r-function can be further weakly reduced to 
satisfy the identity 

dr 2N 



Of 



2fe+l 



= 0, for all k, (4.44) 

{*2i + l=0} 



i.e. tin is requested to be even function of all odd-times £20+1 (this is an example 
of "global characterization" of the weak reduction). Note that (4.44) is imposed 
only on Toda-chain r-function with even values of zero-time. Then (4.44) will 
hold whenever Hi in (4.43) are even (odd) functions of £ dd for even (odd) values 
of i. Since Shur polinomials Pk(t) are even (odd) functions of odd-times for even 
(odd) k, it is enough that the sum in (4.43) goes over even (odd) n when i is 
even (odd). In other words, the restriction on T; m is that 

T lm = fi +m , and f 2k+1 = for all k. (4.45) 

Forced hierarchies. This is another important example of strong reduc- 
tion. It also provides an example of singular r-functions, arising when G = 
ex P {^2 ■Amn'&mrf'n) blows up and normal ordered operators should be used to 
define regularized r-functions. Forced hierarchy appears when G can be rep- 
resented in the form [51] G — GqP + , where projection operator P + is such 



that 

P+ | N) =| N) for N > N , 

(4.46) 

P+ | iV> = for N < N . 
Explicit expression for this operator is 17 

P£=| : exp (- ^ ■= II t 1 " ^) = II (4-47) 

V KN J 1<N KN 

Because of (4.46), P+ | — oo) =), and the identity G | — oo) =| — oo), which 
was essentially used in the derivation in (4.27), can be satisfied only if Go is 
singular and Ti m = oo. In order to avoid this problem one usually introduces 
in the vicinity of sue 



normalized (forced jl r 



all I, m < N 
matrix [52], [51 



i singular points in the universal module space a sort of 
function tL = . One can check that now r/ = oo for 

iv tjv tm 

mill t' can be represented as determinant of a final-dimensional 



t S n = Detff <i,j<N'Hij for N > N ; 



T 



f 
V 



= 1; (4.48) 



r 



I = for N < N Q . 



For N > Nowe have now determinant < if <i //;///< -dinicnsii >nal ( .V - -Am ) \ ( .V - V (l ) 
matrix. The choice of No is not really essential, therefore it is better to put 
No = in order to simplify formulas, phraising and relation with the discrete 
matrix models (No is easily restored if everywhere N is substituted by N — No). 
For forced hierarchies Q can also represent f as 

r f N = Dcto< tJ <N dld{H f , (4.49) 



17 Normal ordering sign : : means that all operators ij> stand to the left of all operators ip. 
The product at the r.h.s. obviously implies both the property (4.46) and projection property 
Pi =P+- 



where — 7ig and d\ = d\ = For forced Toda-chain hierarchy this 
turns into even simpler expression: 

r f N = Bet < i!j<N di+W, (4.50) 

while for the forced Volterra case we get a product of two Toda-chain r-functions 
with twice as small value of N [53] : 



T 2N 



Dct 



o ,., s<t, ' H ' j • (pet < i>j<N di +j (d 2 Hf) 



4[Hf]-4[d 2 Hf]. 



(4.51) 



Forced rj^ can be always represented in the form of a scalar-product matrix 
model. Indeed, 

Kit =^T lm P l _ i {t)P m - j (t) = j j e u ^ +0 ^h l WT(h,h)dhdh, (4.52) 

where T(h, h) = £, m T^h-^h^- 1 , and e u W = e E t >o f ^' = Zi> h l Pl(t). 
Then, since Deto<ij<7vft. 1 = A^(h) - this is where it is essential that the hier- 
archy is forced - 

I ><':,,. vK,, = II/ / e^+^TihMdhidhi ■ A N (h)A N (h), 

(4.53) 

i.e. we obtain a scalar-product model with 

□ 

dfi h - h = e u ^ +0 ^T{h, h)dhdh. (4.54) 

Inverse is also true: partition function of every scalar-product model is a forced 
Toda-lattice r-function. 

KP hierarchy. In this case we just ignore the dependence of r-function on 
times i. Every Toda-lattice r-function can be considered also as KP r-function: 
just operator G KP = Ge H (a point of Grassmannian) becomes i-dependent. 



Usually TV-dependence is also eliminated - this can be considered as a little more 
sophisticated change of G. When N is fixed, extra changes of field-variables 
are allowed, including transformation from Ramond to Neveu-Schwarz sector 
etc. Often KP hierarchy is from the very beginning formulated in terms of 
Neveu-Schwarz (antiperiodic) fermionic fields (associated with principal repre- 
sentations of Kac-Moody algebras), i.e. expansions in the first line of (4.17) are 
in semi-integer powers of z: ipNs( z ) — 'Ylmez ipn,z n ~^ dz 1 / 2 . 

Given a KP r-function one can usually construct a Toda-lattice one with the 
same G, by introducing in appropriate way dependencies on t and N. For this 
purpose t kp should be represented in the form of (4.37): 

r KP {t | G} = Det. iJ<0 r^ p , (4.55) 

where r H.f P = J2iTijPi-i(t). Since T; m is a function of G only, it does not 
change when we built up a Toda-lattice r-function: 

TN{t,t | G} = Detij^Hi+Nj+N] 
Hij =Y J T lm Pi- l (t)P m _ ] (t) = Y / H« p P m ~ J (t). 



(4.56) 



l.rn 



Then 



T KP {t | G} = r {i, | G}. (4.57) 



If we go in the opposite direction, when Toda-lattice r-function is considered as 
KP r-function, 

r {t,t\G} = T KP {t\G(t)}; 
-H« P = J2n> m P m -At) and 

Ti 3 {G(t)} =Y, T i™{G}P m -j(t)- 



KP reduction in its turn has many further weak reductions (KdV and gejissi- 
nesq being the simplest examples). 



4.4 Fermion correlator in Miwa coordinates 

Let us now return to original correlator (4.16) and discuss in a little more 
details the implications of bosonization identity (4.9). In order not to write 
down integrals of J, we introduce scalar field: 18 

(f>(z)= ~^z k +(p + J logz, , 45g . 

fc^O 

kez-o 

such that d(j)(z) = J(z). Then (4.9) states that: 

: VW(A) : = : e^~^ A > : (4.60) 

"Normal ordering" here means nothing more but the requirement to neglect 
all mutual contractions (or correlators) of operators in between : | : \j men Wick 
theorem is applied to evaluate corrletion functions. One can also get rid of the 
normal ordering sign at the l.h.s. of (4.60), then 



= : e*<*> : : e -*< A > : (4.61) 

In distinguished coordinates on a sphere, when the free held propagator is just 
log(z — z), one also has: 

4>{z)i){z) = — ^ : ip{z)ij>(z) : (4.62) 



18 One can consider <j> as introduced for simplicity of notation, but it should be kept in mind 
that the scalar-field representation is in fact more fundamental for generic r-functions, not 
related to the level k = 1 Kac-Moody algebras (this phenomenon is well known in conformal 
filed theory, see [13] for more details). 



Our task now is to express operators e H and e H through the field <j>. This is 
simple: 

H = I U{z)J{z) = <f U(z)d<f>(z) = - <f <f>{z)dU{z). (4.63) 
Jo Jo Jo 

Here as usual U(z) = X^fcx)^ 2 ^ an( ^ integral is around = 0. This is very 
similar to generic linear functional of </>-(A) = — J2k>o I^k^ k , 



H = J 0_(A)/(A)dA, (4.64) 
one should only require that 19 

0U(z) = j> j^jdX, (4.65) 

i.e. 

U(z) = flog (l-j)f(X)d\. (4.66) 
In terms of time-variables this means that 

t k = -^J\- k f(\)d\. (4.67) 

Here we required that U(z = 0) = 0, sometimes it can be more natural to 
introduce also 

|~J, = J log A f(X)d\. (4.68) 

T|hi^ change from the time-variables to "time density" /(A) is known as Miwa 



transformation. In order to establish relation with fcrmionic representation and 
also with matrix models we shall need it in " discretized" form: 



tk — T 
k 



\ a a / 

to = ( E log Xa ~ E log ^° J 

\ a a / 



□ 



(4.69) 



The factor 2iri is included into the definition of contour integral 



We changed integral over A for a discrete sum (i.e. the density function /(A) is 
a combination of S- functions, picked at some points A a , A Q . This is of course 
just another basis in the space of the linear functionals, but the change from 
one basis to another one is highly non-trivial. The thing is, that we selected 
the basis where amplitudes of different (^-functions are the same: parameter £ 



in (4.69) is independent of a. Thus the real parameters are just 



PS 



sitions of the 
ese points in 



points A a , A a , while the amplitude is defined by the density ot-fciic 
the integration (summation) domain. This domain does not need to be a priori 
specified: it can be real line, any other contour or - better - some Riemann 
surface.) Parameter £ is also unnecessary to introduce, because basises with 
different £ are essentially equivalent. We shall soon put it equal to one, but not 
before Miwa transformation will be discussed in a little more detail. 

Our next steps will be as follows. Substitution of (4.64) into (4.69), gives: 

H ^E'MM + ^E'M^)- (4.70) 

a a 

In fact, what we need is not the operator H itself, but the state which is created 
by e H from the vacuum state (N |. Then, since (N | J m = for in < 0, 
(N | e" ?<#, - (A) is essentially equivalent to (N | e~^< A ) with 0_(A) substituted 
by entire </>(A). If £ = 1, e _( ^ A ) can be further changed for ^£>(A) and we obtain 
an expression for the correlator (4.16) where e H is substituted by a product of 
operators ip(X a ). The same is of course true for e H . Then Wick theorem can 
be applied and a new type of determinant formulas arises like, for example, 

r ~ ^t det a ,(JV | ^(A Q )^(A b ) G \ N) (4.71) 

It can be also obtained directly from (4.21), (4.23) and (4.33) by Miwa transfor- 
mation. The rest of this subsection describes this derivation in somewhat more 
details. 



The first task is to substitute 0_ by (f>. For this purpose we introduce operator 

oo 

tkJk = H + + H_, (4.72) 

k— — oo 

where H + = J2k>o^ k ^ k * s J us ^ our °^ H> H- = J2k>o^- k ^ k ^ an< ^ "negative 
times" t-k are defined by "analytical continuation" of the same formulas (4.67) 
and (4.69): 

t-k = \J \ k f{X)d\ = -f ^ - £aA . (4.73) 

Then 

£ t k J k = H + +H-=-dj2<l ) (X a )-J2^a)). (4.74) 

k=— oo \ a a / 

Further j 

e H + +H - = e -i s ^e"+e H - = e^e H ~e H \ (4.75) 
where |— | |^| 

-(*) = E kt ^ = -?Y.l (e - ^ fc ) E ( A * - ^ fe )) = 

fe>0 fc>0 \ a b I 1 / 

2 / , (i-^)(i-^)\ L-l 

|-| (4.76) 
where prime means that the terms wittr a = b are excluded from the product in 
the numerator and accounted for in the infinite "constant", added at the r.h.s. 
In other words, 

e*'M= const JIU^fcMl = 



const • 



n Q iL(Aa-A b ) ; 

f2 (4-77) 

' A 2 (A)A 2 (A) \ 
A(A,A) J 



Since (N | J m = for all m < 0, we have (N | e H ~ — (N \ , and therefore 

(N | e H = (N | e H + = (JV | e H - e H + = e^'® e ff ++ H - . (4.78) 
From eq.(4.74), 



= const • J] : e-^ A «> : : e**< X -> : (4.79) 



where "const" is exactly the same as in (4.77). If £ = 1, eq.(4.61) can be used 
to write: 



{N 1 ^ = A^A^A) (N ' U^WK) (4.80) 



Similarly, 



e n 1 ao = n^nfe) 1 ^T^k' (4 - 8i) 

6 6 A 2 (A)A 2 (A) 



where 



^ = -^E(^-^) 0(4.82) 
and we used the fact that J m | N) = for all m > 0. Finaly, 

TAr{i, t I G} = (iV I e H G e R \ N) 



Hr <„H^n_ A (W A(A,A) 



A2(A)A2(A) A 2(A)A2(A) (4g3) 

■{n 1 nwnw g u^b)i[i>(%) 1 

a a 6 6 

Singularities at the coinciding points are completely eliminated from this ex- 



pression, since poles and zeroes of the correlator are cwweled by those coming 
from the Van-der-Monde determinants. 



Let us now put N — and define normalized r-function 

T „ ( ai G }- T ° (M " |Gl 



(4.84) 



T o{0, I G} ' 

i.e. divide r.h.s. of (4.83) by (0 | G \ 0). Wick theorem now allows to rewrite 
the correlator at the r.h.s. as a determinant of the block matrix: 

/ (0 I V(Aa)V^(A b ) G I 0) (0 I V(A a ) G 4>(% b ) I 0) \ 



dct 



(0 I G I 0) (0 I G I 0) 

(0 I ^(A 6 ) G ip{X a ) I 0) (0 I G V(A Q )<KA b ) | 0) 



V 



(0 I G I 0) 



(0 I G I 0) 



(4.85) 



/ 



Special choices of points A a , . . . , A& can lead to simpler formulas. If A a — > A a , 
so that tk — ► 0, the matrix elements at the right lower block in (4.85) blow up, 
so that the off-diagonal blocks can be neglected. Then 

(0 I e H G I 0) _ 

(4.86) 



T {t,t\G}-,r KP {t\G} = 
A(A,A) 



— det„ 



(0 I G I 0) 
(0 I ^(X a Mh) G I 0) 



A 2 (A)A 2 (A) (0 1 G 1 0) 

This function no longer depends on i-times and is just a KP r-function. 

Matrix element 



¥>(A,A) = 

is singular, when A — > A: ip(\, A) 



(0 I ^(A)V(A) G I 0) 



If now in (4.86) all A — > 00, 



where 



KP 



tp b (X) = (0 I V(A) ( d b -^j ) (00) G I 0) ~ A 



{i|G} 



(0 I G I 0) 
If now fj 

A(A) 



,6-1 



1 + 



(4.87) 



(4.88) 



(4.89) 



This is the main determinant representation of KP r-function in Miwa 
parametrization. 



Starting from representation (4.88) one can restore the corresponding matrix 
H% p in eq.(4.55) [42]: 

= fV-^-V (4.90, 

i.e. 



T« p = iz l ip-i{z). (4.91) 



an 



KP 



Then obviously ^ = H?+ Kj . Now we need to prove that the r-function 
is given at once by dct z y^ Aa - > and Det7i^ p {t}. In order to compare these two 
expressions one should take tk — r Y^a ^a k > so that 



A Q — 2 \ / ^— z — \ a A(A) 

\fe>0 / a=l \ a / a a \ j 

where 

A(A) 



and 



A a (A)= J] (Aa-A g )= ! j 



(4.93) 



< p U E> -* - (n^)E H l(f) Q(A) ^^(A a )- (4.94) 

As far as n is kept finite, determinant of the infinite-size matrix (4.94), 

/KP 



r ::'H,, \ tk= i \- k = ^ SinC(3 ^ 1S 



obvious from (4.94) that the rank of 



the matrix is equal to n. Therefore let usconsidcr the maximal non-vanishing 
determinant, 



n Aa ) det - ( ^VA^f A) )3 e W^(Aa) = (4.95) 



det n jyj(A a ) 
A(A) 



We used here the fact that determinant of a matrix is a product of determinants 
and reversed the signs of i and j. Also used were some simple relations: 



□ 



A a (A) 1 



2 A (A) A2(A)' 
det io jP= (f[X a ) A(l/A), 



A(l/A) = [] f - Y = (-)" ( - 1)/2 A(A) ( [] A. 
thus 



l = l A(A) A* A(A) ' 
Since (4.95) is true for any n, one c an c laim that in the limit n — > oo we recover 



the statement, that T KF {t} = Bet^^jHfj with Hf 3 F given by eq.(4.94) (that 
formula does not refer directly to Miwa parametrization and is defined for any 
t and any j < and i). This relation between (p a 's and Ti.fj P can now be used 
to introduce negative times tk according to the rule (4.58). Especially simple is 
the prescription for zero-time: Tiij — > Tii+Nj+N, when expressed in terms of (p 
just implies that 

det pa(h) det (p a+N (X b ) 



A(A) (detA)*A(A) 
Generalizations of (4.90), like 



(4.97) 



= / / z'z^O | i){z) G i)(z) | Q)e^ k {tkzk+ik? - k) dzdz, (4.98) 



can be also considered. 



4.5 1-Matrix model versus Toda-chain hierarchy 



At the end of this section we use an explicit example of dicrete 1-matrix model 
[54] to demonstrate how a more familiar Lax description of integrable hier^. 
archies arises from determinant formulas. Lax representation appears usual%_ 
after some coordinate system is chosen in the Grassmannian. In the example 
which we are now considering this system is introduced by the use of orthogonal 
polinomials. 

Formalism of orthogonal polinomials was intensively used at the early days of 
the theory of matrix models. It is applicable to scalar-product eigenvalue models 
(see [1] for details about this notion) and allows to further transform (diago- 
nalize) the remaining determinants into products. In variance with both re- 
duction from original iV 2 -fold matrix integrals to the eigenvalue problem, which 
(when possible) reflects a physical phenomenon - decoupling of angular (unitary- 
matrix) degrees of freedom (associated with d-dimensional gauge bosons), - and 
with occurence of determinant formulas which reflects integrability of the model, 
orthogonal polinomials appear more as a technical device. Essentially orthogo- 
nal polinomials are necessary if wants to explicitly separate dependence on the 
the size N of the matrix in the matrix integral ("zero-time") from dependencies 
on all other time-variables and to explicitly construct variables, which satisfy 
Toda-like equations. However, modern description of integrable hierarchies in 
terms of r-functions does not require explicit separation of the zero-time and 
treats it more or less on the equal fooring with all other variables, thus making 
the use of orthogonal polinomials unnecessary. Still this technique remains in 
the arsenal of the matrix model theory 20 and we now briefly explain what it is 



Of course, one can also use this link just with the aim to put the rich and beautifull math- 
ematical theory of orthogonal polinomials into the general context of string theory. Among 

□ 



about. 

In the context of the theory of scalar-product matrix models orthogonal poli- 
nomials naturally arise when one notes that after partition functions appears 
in a simple determinantal form, 

1 n f 

Zn = jvT II / d ^k,h k Bet ikhl ^etjkhi 1 = 

'fc=i (4.99) 

= Dcty J d^ yji^h?- 1 = Detijih?- 1 | h 3 " 1 ), 

(of which eq.(2.53) is a simple example), any linear change of basises h l — > 
Qi(h) = Ylk Aikhk, h? — * Qj(h) = Bjih 1 can be easily performed and 
Z — > Z ■ detA ■ det£>. In particular, if A and B are triangular with units 
at diagonals, their determinants are just unities and Z does not change at all. 
This freedom is, however, enough, to diagonalize the scalar product and choose 
polinomials Qi and Qj so that 



(Qi(h) | Qj(h)) = e+'Sij. (4.100) 

Qi and Qj defined in this way up to normalization are called orthogonal poli- 
nomials. (Note that Q does not need to be a complex conjugate of Q: "bar" 
does not mean complex conjugation.) Because of above restriction on the form 
of matrices A and B these polinomials are normalized so that 

Qi(h) = h l + ...; Q j (h) = h j + ... (4.101) 

i.e. the leading power enters with the unit coefficient. From (4.99) and (4.100) 
it follows that 

i interesting problems here is the matrix-model description of q-orthogonal polinomials. 
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JV 

Z N = Vet 0<lJ < N {V\h?) = JJe*'" 1 . (4.102) 

i=l 

This formula is essentially the main outcome of orthogonal polinomials theory 
fro matrix models: it provides complete separation of the iV-dependence of Z 
(on the size of the matrix) from that on all other parameters (which specify 
the shape of potential, i.e. the measure dfi h this information is encoded in 
a rather complicated fashion in fa. As was already mentioned, any feature of 
matrix model can be examined already at the level of eq.(4.99), which does not 
refer to orthogonal polinomials and thus they are not really relevant for the 
subject. 

Consider now the case of the local measure, dfi h ^ = dfj,hS(h, h), when Qi = 
Qi. The local measure is distinguished by the property that multiplication by 
(any function of) h is Hermitean operator: 



□ 



(hf(h) | g(h)) = (f(h) | hg(h)), if d/x M - ~ S(h - h). (4.103) 
This implies further that the coefficients in the recurrent relation 

i 

hQi(h) = Q t+ i(h) + J2 c vQA h ) ( 4 - 104 ) 

j=o 

are almost all vanishing. Indeed: for j < i 

_ (hQi(h) | Qj{h)) _ {Qi(h) | hQj(h)) _ 



^ {Q^h) | Q {h)) (Q 3 (h) | Qj(h)) 
-^ +1 (Q 3 (h)\ Qj {h))~ ^ 



(4.105) 



In other words, polinomials, orthogonal w.r.to a local measure are obliged to 
satisfy the " 3-term recurrent relation" : 



u 



hQ n (h) - Q n+ i(h) + c n Q n (h) + RnQn-^h) (4.106) 

(the coefficient in front of Q n +i can be of course changed by the change of 
normalization). Parameter c n vanishes if the measure is even (symmetric under 
the change h — * —h), then polinomials are split into two orthogonal subsets: 
even and odd in h. Partition function (4.102) of the one-component model can 



be expressed throng i parameters Ri = of the 3-term relation: 



N-l 

Z N = Z 1 J[ (4.107) 

i=l 

thus defining a one-component matrix model (i.e. particular shape of potential), 
associated with any system of orthogonal polinomials. 

Coming back to the 1-matrix model (2.53), one can say that all the infor- 
mation is contained in the determinant formula (4.37) together with the rule 

V ■ 

H+j- 



(4.34), which defines time-dependence of 7i{j — (h l | h?) = Hi 



dH 



f 

--nL,<=H{^, or 



dt k n *+ k - 



(4.108) 



The possibility to express everything in terms of Tt{ with a single matrix index 
i is the feature of Toda-chain reduction of generic Toda-lattice hierarchy. 

However, in order to reveal the standard Lax representation we need to go into 
somewhat more involved considerations. Namely, we consider representation of 
two operators in the basis of orthogonal polinomials. First, 

n+k , \ uk \ \ 

h k Q n (h) = £ V I 7 = E 1nmQ™{h) (4.109) 

rn=0 y 1 ' m=0 



(here the simplified notation is introduced for (n | f(h) \ to) = (Q n \ f(h) \ Q r , 

, (n I h k I to) 

and 7 W = U 1 M Second, 
(to I to) 



-^7- - - (to i to) QmW = - z^ra^W' 

<tyn = (" I hk I = (fc) 
dt k (n | n) 7ml ' 



(4.110) 



(These last relations arise from differentiation of orthogonality condition (4.100): 
^ d(p n _ d(Q n | Q m ) _ 

dtk dtk n (4.111) 

= I Qm) + (Qn I 8 -%^) + (Qn \ h" \ Q m ) 

at k atk i — i 

by looking at the cases of m < n and m = n respectively.) 

From these relations one immediately derives the Lax-like formula: 

a (fc) n— 1 n+k 



l—m—k l—m+1 



or, in a matrix form, 



where 



9fl = i; -a, 
0*, 



= [i? 7 W, 7 W], (4.113) 



. -Trim if m > n, 
#7^ = <( (4.H4) 
if m < n 



(We remind that usually i?-matrix acts on a function f(h) = J2 n =-oo fnh n 
according to the rule: Rf(h) = J2 n>i f n h n — J2n<i fnh n with some "level" I.) 
These 7 *- fe ^ are not symmetric matrices, but one can also rewrite all the formulas 
above in terms of symmetric ones: 



4» - ei^-^l = ; (4.H5) 

y \m \m)\n \ n) 

From eqs.(4.112) one can easily deduce Toda-equations for 4> n : 

d 2 4> n _ d (n\h l \n) _ 



dt k dti dt k (n | n) 

x in I h k I m) (to I h l I n) 



E-E 



(4.116) 



\m>n m<n/ 



(to I m) (n I n) 

where the r.h.s. can be expressed in terms of R m — e^ m_ In particular, 

Rn +1 - R n = e *-+i-*» - e *»-*»-i. (4.117) 



Let us also mention that in this formalism the Ward identities (Virasoro 
constraints) follow essentially from the relation 

as)' -as (« 18 > 

7 fc>0 

where Hermitean conjugation is w.r.to the scalar product ( | ). For example, 
this relation implies, that 

(Qn | ^) = -(^ | Qn) - E fc ^(Q™ I I Qn)- (4.119) 

fe>0 

Now we note that is a polinomial of degree n — 1, thus (Q n | ^f-) = 0. (In 
fact 

^ = - E f E - - E ^^70 

fe>0 Vm=0 / fc>0 OT ' £ - 1 

Also we recall that (Q n \ h k ^ x \ Q n ) = (Q n | Qn) gt >n l i an( l obtain: 

H fe ^ = (4.120) 

fe>0 
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u 



□ 

for any n. This should be supplemented by relation 5^ = <fi n . In order to get 
the lowest Virasoro constraint (string equation), L-iZn = or L—± logZ/v = 
it is enough just to sum over n from to N — 1. 

For more details about 1-matrix model, Toda-chain hierarchy and application 
of the formalism of orthogonal polinomials in this context see [54]. 

5 r-function as a group-theoretical quantity 

This section contains some remarks about the general notion of r-funcion on 
the lines suggested in ref.[l]. Examples below are taken from [55] and [56]. 

As mentioned in the beginning of the previous section we define the (gener- 
alized) r-function as the generating functional of all the mfljj rix clemeni s of a 
given group element g € G in a given representation R: 



□ 



where operators U and U do no j^-dp pend on R. In order to be even more specific 
one can further request that evolution operators are group elements, i.e. 

AU(t) = U(t) ® U(t) = {U(t) $ I) (I $ U(t)) (5.3) 



r R (t,t\g) = s m,m (*»*)< m\g\m> ( 51 ) 

The choice of functions rll (t, t) is the main ambiguity in the definition of 
r-function and needs to be fixed in some clever way, not yet known in full 
generality. The only a priori requirement is that it is indeed a generating 
functional, i.e. there should be some (^-independent) operators M, acting on 
t, i-variables, which allow to extract all particular matrix elements once T R (t, t) 
is known: 

< m\g\m >=< M R , m ,rh(t,t)\T R (t, t\g) > 

The ambiguity in the choice of i) can be partly fixed (at least in the 

case of the highest weight representations R) by the requirement that 

T R (t,t\g) =< v&c R \U(t)gU{t)\vac R > (5.2) 



where A denotes group comultiplication law. In the case of Lie algebras 
A(Ts) = T$ ® I + I ® T5, and (5.3) is true at least for the evolution oper- 
ators in KP/Toda systems. Later we shall see that in the case of quantum 
groups it can be natural to slightly modify the condition (5.3) |— ] 

Remarkably, the r-function defined in (5.1), always satisfies a family of non- 
linear equations [55], relating t r with different i?'s, which reflect just the fact 
^that matrix elements of the same group element in different representations 
are not independent. Conventional bilinear Hirota equation for KP/Toda r- 
| | functions is nothing but particular case of this generic construction 21 , which 
has two (essentially identical) interpretations: in terms of fundamental repre- 
sentations of GL(oo) and in terms of the level k = 1 Kac-Moody J7(l) algebra. 

5.1 Prom intertwining operators to bilinear equations 

The following construction [55] in terms of intertwining operators is the general 
source of bilinear equations for the r-function (5.1). One can easily recognize the 
standard free-fermion derivation of Hirota equations for KP /Toda r-functions as 
a particular example (with G being the level k = 1 Kac-Moody algebras Gk=i, 
V a fundamental representation, and W - the simplest fundamental representa- 
tion corresponding to the very left root of the Dynkin diagram). Construction 
below involves a lot of arbitrariness. In order to make the consideration more 



21 Note, that it is somewhat different from approach, advocated by V.Kac [49] (see also 
[1]), which makes use of Casimir operators and is less universal than the one to be described 
below (using intertwining operators). 
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transparent, we formulate our construction explicitly for finite-dimensional Lie 
algebras and their g-counterparts. 

Bilinear equations which we are going to derive are relating T-functions (5.1) 
for four different Verma modules V, V, V',V'. Given V, V , every allowed 
choice of V, V' provides a separate set of bilinear identities. Of course, not 
all of these sets are actually independent and can be parametrized by source 
modules V and V and by a weight of finite-dimensional representation. Also 
different choices of positive root systems and their ordering in (5.2) provides 
equations in somewhat different forms. A more invariant description of the 



minimal set of bilinear equations for given G wo dd be clearly interesting to 
find. 



1. Our starting point is embedding of Verma module V into the tensor prod- 
uct V ® W, where W is some irreducible finite-dimensional representation of G 
(in the case of Kac- Moody algebra evaluation representation should be used). 
Once V and W are specified, there is only finite number of choices for V. 

Now we define right vertex operator of the W-type as homomorphism of 
G-modules: 

E R : V — >V®W. (5.4) 

This intertwining operator can be explicitly continued to the whole represen- 
tation once it is constructed for the vacuum (highest- weight) state: 



V=\\n a ) i , = l[(A(T_ a ) n «\0)A, (r,..-„ 

where comultiplication A provides the action of G on the tensor product of 
representations, and 



\{p Q ,j Q } \a>0 



<8 {T- a ) la \Q) v ®\0)w- (5-6) 



For finite-dimensional Ws, this gives every |n a )y in a form of finite sums of 
states |m a )v with coefficients, taking values in elements of W . 

2. The next step is to take another triple, defining a left vertex operator, 

E' L : V' — ► W* ® V, (5.7) 

Note the change of ordering at the r.h.s., this is different from V'<g)W* in the case 
of quantum groups. The product W ® W* of the module W and its conjugate 
contains unit representation of Q. The projection to this unit representation 

tt: W®W* — >I (5.8) 
is explicitly provided by multiplication of any element of W ® W* by 



7T= w(0|® W'{0\1 E 7r {*«'*a} (II ( T +") J °® ( T +«) C J J 



(5.9) 



Using this projection, if it is not occasionally orthogonal to the image of E®E' , 
one can build a new intertwining operator 

r : v®v> E ®S; V ®W®W*®V' I ^ I V®V\ (5.10) 
which possesses the property 

r(s ® g) = [g ® g)T (5.11) 

for any group element g such that 

A(g)=g®g. (5.12) 



1 



3. It now remains to take a matrix element of (5.11) between four states, 
v ,{k'\ v{k\(g®g)r\n)t\n')y, = V ,{W\ v (k\T{g ® g)\n) 9 \n%, (5.13) 
and rewrite this identity in terms of generating functions (5.1). 

5.2 The case of KP/Toda r- functions 

We do not present here the standard derivation of Hirota equations in the free- 
fermion formalism, because it is both well known and easily recognizable in the 
general picture from the previous subsection. Instead we describe here its slight 
variation - starting from the fundamental representations of SX(oo). The rea- 
son why this case is the closest one to the standard integrable hierarchies is that, 
in variance with generic Verma modules for group G ^ SL(2), the fundamental 
representations are generated by subset of the mutually commuting operators, 
not by entire set of generators from maximal nilpoteiQibalgebra. We describe 
the basic construction for G = SL(n), si nce in this case the (finite) Grassman- 
nian construction is the most similar to |thf conventional infinite-dimensional 
(G = 17(1)) situation. 

The Lie algebra SL(n) is generated by operators T±$ and Cartan operators 
Hg, such that \Hg,T±a\ = ±^(a(3)T±s. All elements of all representations arc 
eigenfunctions of Hs, Hs\\) = |(/3A)|A). The highest weight of representation 
F^ k > is jlk- Vectors pk's are "dual" to the simple roots cfj, i = 1, . . . , r: (flidj) = 
S tj , and p = \ J2 S>0 a = fit- 
There are as many as r = rank G = n — 1 fundamental representations of 
SL(n). Let us begin with the simplest fundamental representation F = - 
the n-plet, which consists of the states 

A =T_ (i _ 1) ...r_ 2 r_iV>i, i=l,...,n, (5.14) 



Mor ?o\ er 

T-iipj = Sijipi+i, 

and the weights are given by 



(5.15) 



X(ipi) = pi- ai -...-cU-i, (5.16) 

where pi is the highest weight of F* 1 ). Here T±i = T±£. are generators, asso- 
ciated with the simple roots. Let us denote the corresponding basis in Cartan 
algebra Hi = H 3i , and Hi\X) = |(5iA)|A) = A<|A). Then 

A«=A«(^) = ^(%-^-i)- ( 5 - 17 ) 
This formula, together with (5.14) and defining commutation relations between 
the positive and negative simple-root generators implies that HVMI 2 — L and, 
with the help of the classical comultiplication formula, A(T) = T®7 + 7®T, 
one immediately obeserves that the antisymmetric combinations f/>n . . . ip^i are 
all the highest weight vectors (i.e are annihilated by all A/.(T +i ) and, thus by 
all the Afc(T +( j)). These combinations are the highest vectors of all the other 
fundamental representations F^ , which are thus skew powers of F = fW; 

From this description it is clear that < k < n, moreover and fW are 
respectively the singlet and dual singlet representations. 

According to (5.15) one can also describe all the states of F^ in terms of 
a single generator T_, which is a sum of those for all the r simple roots of G, 

V'i = rr 1 |0)F, i = l,...,n. (5.19) 



Looking at explicit form (5.18) of the states in it is easy to see that they 
can be all generated from the highest weight one, &i2...k by the operators 



R k {Tl) = Ti®I®...®I + I®T l _®...®I + I®I®...®T l _. (5.20) 

These operators obviously commute with each other. For given k exactly k of 
them (with i = l,...,fc) are independent. However, they arc neither (linear 
combinations of) the generators of Lie algebra acting in F 1 -^ , nor even their 
algebraic functions (note that i?fe(Tl) ^ ((i?fe(r_))*). If one wants to make 
clear that F^ k ' is indeed a representation of G, it is better to say that it is 
generated by another set of operators, 



A fe %= T ~s (5.21) 

a:h(a)—i 

The "height" h(a) is the number of items in linear decomposition of the root 
a in simple roors oti (in F^' where all the generators T3 are represented by 
n x n matrices % are matrices with all zero entries except for units at the i-th 
subdiagonal) . In particular, T_ =T\. Operators (5.21) are obviously generators 
of G, instead their mutual commutativit y is somewhat less transparent. Since 
both sets (5.21) and (5.20) generate s a matter of convenience which of 

them is used in particular considerations. In dealing with KP/Toda hierarchies 
the explicitly commuting set (5.20) is more convenient. It is exactly the lack of 
such equivalence of two sets which makes consideration of KP/Toda hierarchies 
more subtle in the quantum (q ^ 1) case, see s.5.4 below. 

The intertwining operators which are of interest for us are 



□ 



F (k) <g> F, 



Here 



r (k) : Fik ~ 1] — " F * ® F {k) , and 
r fc | fe ' : (g> F (k '- 1) — ► F {k) (g> F (k 'K 



F*=F (r) = {^~e« 1 -^ [il ...^ ir] }, 



'(fc) 



^(fe-i) = 9 w ,i 



(5.22) 



(5.23) 



and Tk\ki is constructed with the help of embedding / — > F ® F*, induced by 
the pairing ipi^p 1 '- the basis in linear space ® F^ k induced by Tu^i 

from that in F^ <g> F^'^ is: 



[ti..a k ■ 



□ 



Operation T can be now rewritten in terms of matrix elements 



(5.24) 



as follows: 



(]• 



(k) ( n--^k \ = 
ji ■ ■ ■ jk 



^h...i k \9\^h...j k ) = det _ g 



LT 



,6<fe 



VL71 • - -J&y vjfc+im..-^_i 



□ , 

*x - - - *fc 



<J 



(fc+i) Mi-.-** 



(5.25) 



(5.26) 



Ji ■ ■ • Jk+i 

This is the explicit expression for eq.(5.11) in the case of fundamental represen- 



tations, and it is certainly identically true for any of the form (5. 25). 22 

Let us note that one can use the minors (5.25) to construct local coordinates 
in the Grassmannian. Bilinear Plucker relations satisfied by these coordinates 
are nothing but defining equations of the Grassmannian consisting of all the 
fc-dimcnsional vector subspaces of ?i-dimcnsional vector space. Parametrizing 
determinants (5.25) by time variables (see (5.33)), one gets a set of bilinear 
differential equations on the generating function of these Plucker coordinates, 
which is just a r-function [58]. 

Now let us introduce time- variables and rewrite (5.26) in terms of r-functions. 
We shall denote time variables through Si|r»5|, i = 1, . . . , r in order to emphasize 
their difference from generic tg, is labeled|l55j all the positive roots a of G. Note 
that in order to have a closed system of equations we need to introduce all the 
r times s, for all (though actually depends only on k independent 

combinations of these). 

Since the highest weight of representation F^ k ' is identified as 

\0)fw = l*i fc) fc ^-| [] ( 5 - 27 ) 

we have: 



r {k \t,t\ g) = (* 



i^felexp (X>i4(T[)^ g exp ^>ik(H)) l*Bk>- 

'-(5.28) 



To see this directly it is enough to rewrite the l.h.s. of (5.26) as 



w (ii...i h \ K ( fc '-i) (i'l-'-i'k-A i—i 

(expansion of the determinant gr'* ' in the first column) and now the first two factors can be 
composed into (expansion of the determinant in the first row), thus giving the 

r.h.s. of (5.26). 



Now, 



exp (^URkiT^ = exp (ll k (j2 t i T ^ 



where we used the definition of Schur polynomials 



Essential property of Shur polynomials is that 



Because of (5.29), we can rewrite the r.h.s. of (5.28) as 



□ 



r {k) (t,i\g) 



Px(*)---^(*)<*?i ll2+i2 ,..., fe+i j g 



(5.29) 



(5.30) 



(5.31) 



2+j 2 ,...,k+j, 



■l i 



det ff|(t,t), 

l<a.J3<k ^ 



)P n {t)...P Jk (t) 



(5.32) 



where 



(5.33) 



This formula can be considered as including infinitely many times S{ and Sj, 



and it is only due to the finiteness of matrix g l - £ SL(n) that .ff-matrix is 



I_l 



additionally constrained 



For example, in the simplest case of 



d_\ n 



— I H$ = 0, 



(5.34) 



—H%=0, for i>n. 



The characteristic property of Hp is that it satisfies the following "shift" rela- 
tions (see (5.31)): 



B B 

Tja rra+i rra ttol 



(5.35) 



Coming back to bilinear relation (5.26), it can be easily rewritten in terms 
of iJ-matrix: it is enough to convolute them with Schur polynomials. For the 
sake of convenience let us denote H (^"Jg* ) = dcti< a ,(,<fc . In accordance 
with this notation = H (\' k ^, while bilinear equation turns into: 



[A 



u | ax ...a k ^ H | a k a\ . . . a k _ x j = ^ ^ at x . . . a k [a k ^ ff / <\ , n,, , 



/3i. 



.a fe [a fc 

• flkflk+l 



1 ■ 

ft 



(5.36) 

Just like original (5.26) these are just matrix identities, valid for any Hp. How- 
ever, after the switch from g to H we, first, essentially represented the equationaj 
in the n-independent form and, second, opened the possibility to rewrite them 
in terms of time-derivatives. 



D 



we get: 



fa = i, i = 1, . . . ,k + 1; 
Qij = i, i = 1, . . . , fc — 1, a4 = fc + 1 
$ = i, i=l, ...,fc-l 

□ 

H ( 1 --- k \ H (k + l,l...k-l 
l...k) \k + l,l...k-l 



l...k-l,k \ (k+l,l...k-l 
l...fc-l,fc + l 



□ 



(5.37) 



fc, l...,fc-l 



(5.38) 



(!: 


.AH 


-1^ 


i-C: 


.fc- 


!) 




.AH 


-1/ 




.fc- 





(all other terms arising in the process of symmetrization vanish). This in turn 
can be represented through r-functions: 

dl d lT (k) . r « _ d lT W dT {k) = r (fc+i) T (fc-i). (5.39) 
This is the usual lowest Toda-lattice equation. For finite n the set of solutions is 



labeled by g S SL(;i) (rather than SX(oo)) as a result of additional constraints 
(5.34). 

We can now use the chance to illustrate the ambiguity of definition of r- 
function, or, to put it differently, that in the choice of time- variables. Eq.(5.39) 
is actually a corollary of two statements: the basiq-4dentity (5.26) and the 



particular choice of evolution operators in eq.(5.2), wfeidh in the case of (5.28) 
implies (5.33) with P's being ordinary Schur polynomials (5.30). At least, in 
this simple situation (of fundamental representations of SL(n)) one could define 
r-function not by eq.(5.2), but just by eq.(5.32), with 
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H<$ (t, t) — > ?^ (t, t) = E ^ ^ © (5.40) 



with any set of independent functions (not even polynomials) V a - Such 



_(*0 



det H% (5.41) 



still remains a generating function for all matrix elements of G = SL(n) in 
representation F^ k ' . This freedom should be kept in mind when dealing with 
"generalized r- functions" . As a simple example, one can take V a {t) to be q- 
Schur polynomials, 

Y[e q {t^) = Y J Pf{t)z\ or 

(5.42) 



which satisfy 



D u Pf{t) = {D ti yp) q \t) = P) q \{t). (5.43) 
where D are finite-difference operators. Then instead of (5.35) we would have: 

D ti H1 = H% + \ D- ti H%=Hl +i (5.44) 

and 

r^(t,t\g) = det D?- 1 D^ X U\ (t,t). (5.45) 

l<a,p</c 1 

So defined r-function satisfies difference rather than differential equations [59, 
60]: 



□ 



□ 



□ 

r« • D tl D- tl r^ - D tl r<*> ■ D h r^ = r^- 1 ) • M+Mj+T^+V, 

(5.46) 

We emphasize, however, that this is just another description of the SL(n), 
not SLq(n) r-function, if it is interpreted as a generating function of matrix 
elements. In particular, this r-function takes c- rather than g-number values. 
Still, as concerns its times-, not g-dependence, it has something to do with the 
SL q (n) group, in the spirit of relation between g-hypergeomct jric| functions and 
quantum groups (see, for example, [61]). 



5.3 Example of SL{2) q 




Construction from subsection 5.1 is immediately applicable to the case of quan- 
tum groups, the only thing one should keep in mind is that our definition (5.1) 
gives r as an element of "coordinate ring" A(G q ), not just a c-number. If one 
wants to obtain a c-number r-function for q =/= 1 it is necessary to restrict the 
construction further to particular representation of coordinate ring (this last 
step will not be discussed in this paper). We present here in full detail the 
simplest possible example of SL(2) q [55]. 
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5.3.1 Bilinear identities 

To begin with, fix the notations. We consider generators T+, T_ and To of 
U q (SL(2)) with commutation relations 

□ 



□ 

q To T±q -n = 5 ±i T±) 



q2T _ q—2T 
q-q- 1 



(5.47) 



and comultiplication 



(5.48) 



A(T±) = g T ° <g> T± + T± <8> g _T ° , 
A(g T °) = q To (g) g T °. 

Verma module V\ with highest weight A (not obligatory half- integer), consists 



of the elements 



such that 



|n) A = r l |0) A , n>0, 

T_\n) x = \n + l) x , 
To\n) x = (A - n)\n) x , 
T + \n)\ = b n (\)\n - 



(5.49) 



(5.50) 



6 n (A) = [n][2A+l-n], [x] = 



IMIa = *( n \ n )x 



q - q~* 

[n}\ T q {2\ + 1) agz/2 [2A]![n]! 
T q {2X + l-n) 



[2A- 



Now, 



(A(T_)) n = q nTo (g) T" + [n]T_g ( "- 1)To <g) T^q-^ + ...+ 
+ [n]TT 1 q To ® T-q- {n -^ T « + T™ <g> g-" T °. 



(5.51) 



Let us manifestly derive equations (5.13) taking for 14 7 an irreducible spin-i 
representation of U q (SL(2)). Then V" = V^i, V ~ Vx and the highest weights 



of V" in W ® V or 1/ <g> are 23 : 

|0) A+J =|+)|0) A) \+) = \o) i , 
or |0) A |+); 

|0>a- 4 = |+)|1)a -9 (A+i) [2A]|-)|0) A , I-) = |l)i, □ 

or {q^q- 1 ) | 1) x |+) - g"< A+ *> [2A] |0) a R . 
Entire Verma module is generated by the action of A(T_): 

|n) x+i = (A(T_)) n |0) A+ i — » 
q n/2 (\+)\n)x+q- {X+ ^[n]\-)\n-l)x), 
or g -»/a(|„) x |+) + g(*+*>[ n ]| n _l) x |_; 



(5.52) 



(5.53) 



(5.54) 



\n)x_ h = (A(T_)) n |0> A _ 4 — > 
g "/ 2 (|+)|n+l)A + g (A+ ^[n-2A]|-)|n)A) , (5.55) 
<T" /2 (|n + 1)a|+> + g- (A+ ^[n - 2A]|n) A |->) ; 



Step 2 to be made in accordance with our general procedure is 
tensor product of two different W's onto singlet state S = |+)|— ) — q 

(A\+) + B\-)) ® (\+)C + \-)D) — ► AD — qBC. 
With our choice of W we can now consider two different cases: 



to project the 

i \ .24 



(5.56) 



'Hereafter we omit the symbol of tensor product from the notations of the states |+) (g> |0); 



etc. 



24 This state is a singlet of U q (SL(2)). In the case of U q (GL(2)) one should account for 
the U(l) non-invariance of S. This is the origin of the factor det q g at the r.h.s. of the final 
equation (5.64). 
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(A) both V = V x _i and V' = V x ,_i, or 

(B) V=V x _ l2 and V' = ^, + i: 
Case A: 



\n) x _i\n') x ,_i — > 
([n'-2A'] g A '|n+l) A |n')v ~ [n - 2\]q- x \n) x \n' + l) x , 



(5.57) 



Case B: 



\n) x+ i\n') x ,_i 



q^ 1 ([ri - 2X']q y \n) x \n' ) v - [n]g +A+1 |n - l) A |n' + 1) A ,) 



(5.58) 



Now we proceed to the step 3. Consider any "group element" , i.e. an element 
g from some extension of U q (G), which possesses the property: 



A{g)=g(g>g, 
and take matrix elements of the formula (5.11): 

\>{k'\ x(k\(g®gT = T g®g)\n) x \n') x , 
The action of operator T can be represented as: 

i,v 

and in these terms (5.60) turns into: 



(5.59) 



(5.60) 



(5.61) 



T(k,k'\m,m') 



" m ir(m\g\n) x ( m '\g\n') xl 



\m - \\m ,, 
i ii A n n A 



(5.62) 



= J2{k\9\l)x{k'\g\l')x'T(l,l'\n,n'). 



Li> 



In order to rewrite this as a difference equation, we use our definition of 
T-function: 



T X (t,t\g) = (X\e q (tT+)ge q (iT-)\X) = ^(m| 5 |n) A 



TTlrTTT" 



(5.63) 



m ! \n\ 

Then, one can write down the generating formula for the equation (5.62), using 
the manifest form (5.57)-(5.58) of matrix elements T(l, l'\n, n'): 

Case A: 



MfM+ ( q y Df ) t'D^' ) - q-HDf X) DP) T X (t,i\g)T X ,(t',t'\g) = 

= [2\}[2\'}(det q g)^M^M+ («r<*+*V - q^'+^t) r x _ h {t,t\g)r xl _ h {t' ,¥\g). 

(5.64) 

Here d[°^ = - — ^ r ij t M ' an^~~|lf ± are multiplicative shift operators, 
Mff(t) = f(q±H). 



Case B: 



jMfM+ (q^'t'D^ - q^HDf) T X (t^g)ry(t' ,¥\g) 

[2A ' ] M> M+ (q* tD? x ^ - qH> D^) T X+i (t^g)T x ,_yti,¥\g). 

(5.65) 



[2A + 1] 



Let us note that the derivation of these equations can be presented in the form 
which looks even closer to conventional frcc-fcrmion formalism. It is possible to 
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represent operator T in component form as (g) E^ — qE^ <S> E^, where Ei's 
are components of the vertex operator (given by fixing different vectors from 
W). Then the equation (5.11) can be rewritten 

^<0|e g (tT+) E?g e g (tT_)|0)y -y, (0\e q (tT+) E^g e q (iT-)\0)v> - 
-q v (0\e q (tT+) E?g e q (tT_)\0) v - v , (0\e q (tT + ) E{g e q (tT_)\0) v , = 

(5.66) 

= y(0\e q (tT + ) gE* e q (tT_)\0) v - v , (0\e q (tT+) gE^ e q {iT_)\G) v , - 
-q v (Q\e q (tT + ) gE* e q (iT_)\0) v - yl (0\e q (tT + ) gE{ e q (iT_)\0) v , 

We can easily obtain commutation relations of E^s with generators of algebra 
as well as their action on vacuum states. Then, it is straightforward to commute 
E^s with g-exponentials in the expression (5.66) and represent the result of the 
commutation by the action of difference operators. Of course, the results (5.64) 
and (5.65) are reproduced in this way. 

5.3.2 Solution to bilinear identities 

In this particular case (of SL(2) q ) one can easily evaluate the r-function explic- 
itly, and let us use this possibility to show how bilinear equations are satisfied. 
The fact that our r-function is operator-valued jvil^ be of course of principal 
importance. 

Let us begin from the case of A = i. Then 



ri(t,t\g) = (+M+) +i(+\9\-) + t(-\g\+) +tt(-\g\-) = 
= a + bi+ct + dti, 



(5.67) 



where a, 6, c, d are elements of the matrix 

/ a b 

T = | I (5.68) 

with the commutation relations dictated by TT1Z = 1ZTT equation [62] 

□n 

ab = qba, 



□ 



□ 



ac = qca, 
bd = qdb, 
cd = qdc, 
be = cb, 



(5.69) 



ad — da = (q — q~ 1 )bc. 
If 6 or c or both are non-vanishing, Ti(t,t\g) with different values of time- 
variables t, i do not commute. Still such ri {t, t\g) does satisfy the same bilinear 
identity (5.64), moreover, for this to be true it is essential that commutation 
relationsrfe.69) are exactly what they are. Indeed, the l.h.s. of the equation 



(5.64) is 



equal to 



-qi^Mf {b + dt)^j~M^{a + ct') + q~i ^Mf (a + ci)y / M? (b + dt') = 
= (q~^ab - qHa) + (q~?cd - q? dc)tt' + (q~^cb - q^ da)t + (q~^ad - q~Hc)t' 

= {q-^t' -qh)det q g, 

(5.70) 

which coincides with the r.h.s. of the equation (5.64). 

To perform the similar check for any half-integer-spin representation, let us 
icfeponding r-function can be easily written in terms of ri . 



note that the co 



Indeed, 



\n)x = {q To (8 T_ + T_ (g) <r T °) n |0) A _i <g> |0)i = 



n-1 



\\ n \ = A-± 

-n/2 



= (q T ° <g) T+ + T + 



(5.71) 



?-" /2 (x-i<»l® i(0| + [% A A _ 4 (n-l|® 4 (1|). 



□ 



Thus 



□ 



□ 

Using the relations like 

q?ab + q~Ha = [2}qHa = [2]q~iab etc. (5.76) 
one gets for this case 

Ti(t,t\g) = a 2 + [2}q-?abi + [2]q^^ act + b 2 P + {[2]qbc + [2}da)ti + c 2 t 2 + 
+[2]qidbtf + [2}q?dct 2 i + dH 2 ?. 

(5.77) 

With this explicit expression, one can trivially make the calculations similar to 
(5.70) in order to check manifestly equation (5.64) for A = 1, A' = i, 1 and 



+ q x [k] A-i(fe - 1|<?K>a-4<-|<?I+) + r A WN x-l( k - Mgln- lK-lHsR equation (5.65) for A = A 



(5.72) 

or, in terms of generating (r-) functions: 

T X (t, t\g) = y/MfMf (r A _j (t, i\g) (a + q X tb + q X tc + q 2X tid)) . (5.73) 
Applying this procedure recursively we get: 

T X (t,i\g) =T x _i(q-h,q-n\g)Ti(q X -h,q X ~n\g) = 



M _ 1 



if AeZ/2 



Ti(qi- X t, q^- X t\g)Ti{qi-% q^ x t\g) ...Ti {q x ~H, q x -?i\g) 



A-* 



(5.74) 



i.e. for half-integer A r A is a polynomial of degree 2A in a, b, c, d. 
For example, 

n(t,t\g) = Ti(q-*t,q-$i\g)Ti(qh,q$i\g) = 
— (a + q~Hb + q~^tc + q~ X tid)(a + q^ib + q^ic + qttd) = 

— a 2 + (q^ab+ q~Ha)i + (q? ac + q~ ' ca)t + b 2 P+ 
+(qad + bc + cb + q~ 1 da)ti + c 2 t 2 + (qHd + q^^db)tt 2 + 
+ (qicd + q~idc)t 2 i+ d 2 t 2 P 



(5.75) 



Thus, we showed explicitly (for the case of SL q {2)) that the quantum bilinear 
identities have as many solutions as the classical ones, provided the r-function 
is allowed to take values in non-commutative ring A(G). 

5.4 Comments on the quantum deformation of KP/Toda 
r-functions 

As we saw in the previous subsection, the generic construction is easily appli- 
cable to quantum groups, still the problem of quantum KP/Toda hierarchies 
deserves separate consideration and is not yet fully resolved. The problem is, 
that the evolution operator U{t) in (5.2) is usually constructed from all the oper- 
ators of the algebra, not just from a commuting set - as it happens in particular 
case of fundamental representations of GL(N). As result generic evolution of 
r with variation of i's is not described as a set of commuting flows, rather they 
form a closed, but non-trivial, algebra. This manifests itself also in the fact 
that naturally the number of independent time-variables is rather close to di- 
mension than to the rank of the group. The problem of quantum deformation of 
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KP /Toda hierarchy is to find a deformation which, while dealing with r-function 
for quantum group, is still describable in terms of few time- variables. If at all 
resolvable this is the problem of a clever choice of the weight functions s^ m (t, t) 
in (5.1). Following [56] we shall now demonstrate that the problem is resolvable 
in principle, though at the moment it is a quantum deformation of somewhat 
non-conventional description of KP/Toda system (with evolution, introduced 
differently from that in s.5.2, and it is not just a change of time- variables: the 
transformation is representation i?-dependent). 

According to [57] parametrization of group elements which allows the most 
straightforward quantum deformation involves only simple roots ±c?j, i = 
1, . . . ,r G : 

9 = 9U9D9L, 

r G (5.78) 
9u =l[e^, 5L =n eXsT -* <3) ' 9D = i[e m 

3 S i—1 

Every particular simple root cfj can appear several times in the product, and 
there are different parametrizations of the group elements of such type, depend- 
ing on the choice of the set {s} and the mapping i(s). Quantum deformation 
of such formula is especialy simple because comultiplication rule is especialy 
simple for generators, associated with simple roots: 

A(T 4 ) = T t ®q- 2H > +J®T i , 

□ (5.79) 

A(T_j) = T_, ® I + q 2H > ® T_i 

□ 

For q 1 any expression of the form (5.78) remains just the same, proviFfeJ 
exponents in gjj and g^ are understood as g-exponents (in the simply-laced case, 
/ 2 -exponents in general), and parameters "0,x,</> become non-commuting 
generators of "coordinate ring" A(G q ). Actualy thepr-form a kind of a very 



simple Heisenberg-like algebra: 



XsXs> = q S * <3)S ^'>Xs'Xs, s < s', 

(5.80) 



= e s q M q^\ 



q^Xs = Xsq^q^ 

These relations imply that A(<?) = g ® g. 

The simplest possible assumption about evolution operators would be to say 
that, just as it was in the case of the standard KP/Toda theory (see s.5.2), U(t) 
is always an object of the type gu, while U(i) - of the type g^. However, these 
are no longer group elements: 

Mgu)^9u®gu, A(g L )^g L ®g L , 

because of the lack of factors go- Still the simplest possibility is to insist on 
identification of U and U as objects of the type gu and gu respectively, and 
explicitly investigate implications of the failure of (5.3). As result one obtains 
instead of (5.3) 



□ 



A(u(0) = uP(0-u^(0,^ ( 5 - 81 ) 



where 



u(o= n < 

s 

= II < £ * & T *« ® t~ 2H%{s) ) *I®U®, ( 5. 8 2) 

S 

u ( b ] = n < £ i i^ 1 ® t ^)) =i® u(o 

s 

and this has some simply accountable implications for determinant formulas for 
quantum r-functions. 



In s.5.2 we essentially used an evolution operator of the type 

r 

U ® = IT IT expfeT^) (5.83) 

l<i<AT l<i<i 

where ^ are certain functions of only TV idependent variables t. While (5.83) is 
trivial to deform in the direction of q ^ 1, it is a separate (yet unresolved) prob- 
lem to find such reduction to only iV- variables, consistent with the commutation 
relations between 6j , 

iiiii'f = q'^-^'-i'&j'Zij, {hj} < {i'J'}- 

One can instead use a much simpler evolution, 

TO < 

u(o = n ex p te^) ( 5 - 84 ) 

i=l 

This is enough to generate all the states of any fundamental representation from 
the corresponding vacuum (highest vector) state, but < vac Fn \ U{^) has nothing 
to do with the usual < vac Fn \ U(t), where U(t) is given by (5.83). It can be 
better to say, that identification < vac Fn \ U^ A \^) = < vac Fn \ U{t) defines a 
relation £i(i), which explicitly depends on n. 

One can of course built the theory of KP/Toda hierarchies in terms of £- 
variables instead of conventional i-variables, but it can not be obtained by just 
change of time- variables: the whole construction will look different. Instead 
this new construction is immediately deformed to the case of q ^ 1: instead of 
(5.84) we just write 



i=l 

where £'s are non-commuting variables, 



66 </ " " %<j, (5-86) 



and-it is easy to derive quantum counterpart of any statement of the classical 



(q^^l) theory once it is formulated for £-parametrization. 

In what follows we first briefly describe the conventional KP /Toda hierarchy 
in this non-standard parametrization, then consider the corresponding quantum 
deformation and derive the substitute of determinant formulas for r„ = T Fn in 
the case of q ^ 1. 



,0 On 



the modified KP/Toda hierarchy 



Our first purpose is to demonstrate that all the main ingredients of description of 
the classical KP/Toda hierarchy, as described in s.5.2, are preserved if evolution 
(5.85) is used instead of (5.83), in particular, there are determinant formulas 
and a hierarchy ol differential equations. 



From now on we denote the r-function associated with the evolution (5.85) 
through T(£,£|g). This r-function is linear in each time-variable £i, hence, 
it satisfies simpler determinant formulas and simpler hierarchy of equations. 
Indeed, now we have 



h(t,Z\g) = < o Fl \U(Z)gU(d\o Fl > = Sk7s * < *M* > (5.87) 

k,k>0 

where s k = so = 1, and 



fr™(£,£|<?) = < m Fl \U^)gO(0\fh Fl > = - V s k s- k < k\g\k > 



k>f. 



^ ^ Q Q _ 1 d d - 

— — ^ T\ «1 — T i(£>£lff) = = -^-^ Tl (^£\9)- 

k>m 

(5.88) 

□ 



Thus, 25 



tn+i = det t[" 

< m . rn < n 



n 



dot,-, 



s n s 



k,k> 7i 



det 

0<m,m<n— 1 



V 9km 9mm 



sfeSfc < k\g\k > 



k2 



s n s 



/ 

(n) 



k,k~>7 



(5.89) 



5.4.2 q-Determinant-like representation 

In this section we demonstrate how the technique developed in the previous sec- 
tions is deformed to the quantum case and, in particular, obtain q-determinant- 
like deformation of (5.89). Our evolution operator (5.85) satisfies the following 
comultiplication rule: 

n 

A"" 1 (U{T t }) = Y[ (5.90) 

m— 1 

where 

£/(™) = u { I ® I <8> 6Ti ® q~ 2H > ® . . . (8 q~ 2ffi } (5.91) 

25 One can compare determinant representations (5.41) and (5.89) to find the connection 
between different coordinates t and £. For every given n the variables are some functions of 
Pj(t). For example, in the simplest case of the first fundamental representation i^ 1 ) we have 
T l(*l9) = ft-(Cls) an d s/c = Pk(t}> Jjf - = X/i s '— k '&7~ s However, identification of r n (i) and 
Tn(£) with n^l will lead to different relations between f; and t. Thus the two different evolu- 
tions are not related just by a change of time-variables, relation is representation-dependent, 
and can not be lifted to the actual KP/Toda c ^se | (when n = oo). Two evolutions provide 
two equally nice, but not just equivalent descriptions of the same hierarchy. 



□ 

and Ti appears at the m-th place in the tensor product. Similarlyj 

t/( m ) = {J { q 2H * ® q 2H i g> T_ 4 ® I® . . . ® 1} . (5.92) 

Now let us transform the operator-valued q-factors into c-number ones. Let 

Hi\jF 1 >= \j\3Fi >, < 3Fi \Hi = hi,j < j Fl | 
(in fact for SL(N) = +1, 2hi t i = — 1, all the rest are vanishing). Then 

= (®IU < Jm|) & n -\U)g® n /\ n -\U) (SjUlim >) = 

71 

II < Jm\U {^T iq - 2 ^'=^ } g U {^T.^ZT } \] m > = 



m= 1 



m=l 

1 (5.93) 



In order to get the analogue of (5.41), on b sh ould replace antisymmetrization 
by q-antisymmetrization, since, in quantum case, fundamental representations 
are described by q-antisymmetrized vectors. We define q-antisymmetrization as 
a sum over all permutations, 

([!,..., k] q ) = 5>q) dog P (P(l), . . . , P(k)) , (5.94) 

□ 

where 

deg P = # of inversions in P. (5.95) 
Then, q-antisymmetrizing (5.93) with jk = k — 1, j% — k — 1, one finally gets 



deg P+deg P' . 



X 
P,P' 

(5.96) 

m=0 ^ ' 

This would be just a (/-determinant, be there no g-factors which twist the time- 
variables. 26 

To make this expression more transparent let us consider the simplest example 
of the second fundamental representation: 

T2 = T™(qZ 1 ,q- 1 &,fr, 6,6,6)^(6,6,6 <?6,6,6)- 
-rtHtti, <T X 6,6; 6,6, 6)n 10 (6,6, 6 T^i, <z6,6)- 

(5.97) 

6,6,6)^(6,6,6 <?6,6,6)+ 
+9M 1 (9" 1 6, 6, 6 6,6,6K°(6,6, 6; 9 _1 6, ?6, 6)- 

This can be written in a more compact form with the help of operators 

H>f = Di ® I, Df = M7 SiS * ® Dj, 

of = A ® II A V 3l3j . ®f = 1 ® A. 

These operators have simple commutation relations: 

Bfof = q s ^Bf Il)f . 



(5. 



(5.99) 



26 Let us remind that the g-determinant is defined as 

dct 9 A ~ Al . . . = ^ (-^ ^ P' [J A JJJ) 

P,P' a 

Note that this is not obligatory the same as . . . A"j . It is the same oCy-lfor peculis 
commutation relations of the matrix elements A 3 - . 



Then, 



r 2 = (Mf M+) (Df - gDf ) • (m>f - gD^ n <X> t x . (5.100) 

6 Conclusion 

These notes combine presentation of some well established facts with that of 
more recent and sometime disputable speculations. There are all reasons to 
believe that further developements will prove that the theory of generalized 
r-functions and non-perturbative partition functions can become a flourishing 
branch of mathematical physics with applications well beyond the present mod- 
est scope of topological theories and c < 1 string models and with profound 
relations to other fields of the string theory. It is also important that there are 
plenty of "small problems" at all the levels of this theory, which are enjoyable 
to think about. 
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